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Abstract 



The adiabatic particle number in mean field theory obeys a quantum Vlasov 
equation which is nonlocal in time. For weak, slowly varying electric fields this 
particle number can be identified with the single particle distribution function in 
phase space, and its time rate of change is the appropriate effective source term 
for the Boltzmann- Vlasov equation. By analyzing the evolution of the particle 
number we exhibit the time structure of the particle creation process in a constant 
electric field, and derive the local form of the source term due to pair creation. 
In order to capture the secular Schwinger creation rate, the source term requires 
an asymptotic expansion which is uniform in time, and whose longitudinal mo- 
mentum dependence can be approximated by a delta function only on time scales 



much longer than \Jp\ + m 2 c 2 /eE. The local Vlasov source term amounts to 

a kind of Markov limit of field theory, where information about quantum phase 
correlations in the created pairs is ignored and a reversible Hamiltonian evolution 
is replaced by an irreversible kinetic one. This replacement has a precise counter- 
part in the density matrix description, where it corresponds to disregarding the 
rapidly varying off-diagonal terms in the adiabatic number basis and treating the 
more slowly varying diagonal elements as the probabilities of creating pairs in a 
stochastic process. A numerical comparison between the quantum and local ki- 
netic approaches to the dynamical backreaction problem shows remarkably good 
agreement, even in quite strong electric fields, eE ~ m 2 c s /h, over a large range 
of times. 



* Deceased 



I. INTRODUCTION 



In recent years there has been considerable interest in establishing the precise connection 
between quantum field theory and classical kinetic theory. This interest is motivated by the wide 
variety of problems in different fields of physics which require a consistent description of quantum 
many body phenomena far from equilibrium. Examples include chiral symmetry restoration and 
the quark-gluon plasma phase of QCD, soon to be probed by relativistic heavy-ion colliders, 
baryogenesis at the electroweak phase transition, and the formation and decay of topological 
defects or Bose condensates, whether in the hot, dense early universe, or a cryogenic laboratory 
environment. 

At their root all these systems may be treated as field theories with well-defined Hamiltonian 
evolutions and (except for the case of explicit CP violation in the electroweak theory) micro- 
scopic time reversal invariance. Yet, a large body of experience confirms the macroscopically 
irreversible behavior of such systems far from equilibrium, so that it should be possible to approx- 
imate the unitary Hamiltonian evolution of such systems by an irreversible kinetic description, 
under suitable circumstances. In addition to the numerous potential applications, this raises the 
fundamental issue of the precise connection between microscopic reversibility and macroscopic 
irreversibility which lies at the heart of much of nonequilibrium statistical mechanics. 

Whereas kinetic theory is by far the more developed and familiar framework to study nonequi- 
librium problems, the first steps in their practical solution in the context of quantum field theory 
have been taken only relatively recently Until that had been done it was difficult to even 

formulate the question of the relationship between the field theory and kinetic theory approaches 
in a clear way. As a practical matter the kinetic description is certainly the simpler one to 
formulate and implement numerically on a computer. However, the Boltzmann-Vlasov equation 
essentially describes classical point particles, and extensions to quantum collective phenomena, 
time-evolving mean fields and off-shell virtual processes, which are quite natural in field theory, 
present considerable difficulties for a purely kinetic approach. Also lost in the kinetic description 
from the very outset is a detailed understanding of how time reversible Hamiltonian evolution 
comes to be replaced by time irreversible dissipative behavior. For these reasons of both fun- 
damental interest and practical application, our purpose in this paper is to explore the precise 
relationship between the two approaches in a concrete example. 

In the interest of being as clear and specific as possible we focus our attention in this paper 
on charged particle creation in electric fields, a phenomenon which was discussed nearly seventy 
years ago by Klein and Sauter, and twenty years later by Schwinger as a prime example of the 
then newly developed theory of quantum electrodynamics || . Over the years there has developed 
an extensive literature on the topic which has continued to attract interest up to the present 
time [P|-|T7|]. Several monographs summarizing this activity have also appeared \TE\. Given this 
background it might be supposed that no aspect of particle creation in electric fields has been 
left uncovered. The reason that this is not quite true is that interest in the real time evolution 
of particle creation and its incorporation into transport theory by an effective source term is 
relatively recent. 

In many treatments of particle creation, analytic continuation of the amplitudes to complex 
time have been employed. Though elegant and useful in other contexts, complex continuation 
methods cannot address directly the real time evolution of the particle creation event and thus 
cast little light on the source term for a kinetic description. A suggestion of how to incorporate 
the Schwinger pair creation mechanism in the context of kinetic theory was first made in 1979, 
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based on an intuitively appealing picture of the instantaneous semiclassical creation event [ID]. 
This mechanism has been a subject of renewed interest in the context of heavy ion collisions and 
QCD due to the suggestion that the receding ions might produce a strong chromoelectric flux 
tube between them which shorts itself out by the creation quark/ ant i- quark pairs (see e.g. @ 
and references therein). The ansatz of ref. ]IDJ has been taken over to the QCD flux tube model 
as well. Yet it should be clear from the outset that a delta function source term which requires 
that the charged particles be created at precisely zero momentum, at a definite instant of time 
can only be an approximation to the rapid but continuous evolution of wave amplitudes in the 
underlying quantum theory. Calculations of the backreaction of the charged particle pairs on the 
electric field in QED in a well defined continuous evolution were compared with the ad hoc kinetic 



theory, according to the ansatz of ref. [ lOj . Reasonable qualitative agreement between the mean 



field evolutions in the two approaches was found, although they certainly differ in quantitative 
detail, such as in the distributions of created particles @. In these numerical investigations the 
time structure of the individual creation events was not addressed, leaving open the question of 
the limit of validity of the delta function ansatz for the source term. 

The Wigner function formalism has also been proposed [I^J2pj as a method for deriving 
relativistic transport equations from the underlying field theory. It has become increasingly clear 
however, that the covariant Wigner function does not readily lend itself to practical calculation, 
because covariance requires splitting the time variable in the Wigner transformation in parallel 
to the splitting of the spatial variable, with the consequence that the problem ceases to be well 
posed as an evolution from initial data. More recently, an alternative, noncovariant formalism, in 
which the time variable is not split has been suggested |^3[l4[ . As has been emphasized in earlier 



work the lack of manifest covariance is not a problem since the initial value description of even a 
relativistic field theory in Hamiltonian terms is necessarily noncovariant in form, but the evolution 
equations are completely equivalent to those derived from a covariant action principle @ . In any 
case, a firm conclusion about the source term has not been obtained by these investigations 
either. Finally, the general projection formalism of Zwanzig |21] has been advocated as a route 
to a transport description of particle creation [T5|-|iTf , although the time structure of the creation 
process itself has not been investigated in detail in this approach, and the conditions of validity 
of the delta function approximation for a local source term in the Vlasov equation has remained 
obscure. 

By revisiting the electrodynamic pair creation problem our purpose in this paper is to elu- 
cidate fully the precise connection between the field theoretic and kinetic treatments in this 
particular case. Application and extension of our methods of incorporating particle creation into 
a kinetic description for the other situations of interest will then become possible. Our first step 
will be to specify completely the adiabatic particle number basis in which particle creation can 
be described as a phase interference (or dephasing) phenomenon of the quantum theory from 
the effective Hamiltonian point of view ||. Writing the explicit Bogoliubov transformation to 
this adiabatic particle basis then identifies a time dependent particle number whose total change 
recaptures the Schwinger formula in a constant, uniform electric field, and whose time derivative 
yields the appropriate source term for the Boltzmann- Vlasov equation. The adiabatic particle 
number obeys a nonlocal quantum Vlasov equation, and in this sense is completely consistent 
with the general approach advocated in refs. p^j— [T7]] . The relationship of our method to that 
of the projection formalism may be seen most clearly by considering the density matrix in the 
adiabatic particle number basis. However we have no need for the general projection formalism, 
since the source term for the Vlasov equation can be written in closed form in terms of the wave 
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functions of the charged particle modes in the background constant electric field. In this way we 
derive for the first time a local form for the source term, which explicitly exhibits the relationship 
to the semiclassical picture of particles spontaneously appearing out of the vacuum in real time. 
The electromagnetic current of the charged particle pairs also has a simple form in this basis, 
corresponding to a clear physical interpretation in terms of a quasiclassical conduction current 
and the quantum polarization current of particle creation. The fact that the current grows lin- 
early in time for a fixed external electric field and that therefore backreaction must eventually 
become important even for arbitrarily small coupling is also easy to see in the adiabatic particle 
basis. This will also serve to clarify the nature of the "time divergences" discussed in refs. [f|. 

The essential physical ingredient in passing from the quantum unitary evolution to the irre- 
versible Vlasov description is the dephasing phenomenon, i.e. the near exact cancellation of the 
rapidly varying phases of the quantum mode functions contributing to the mean electric current 
of the created pairs. This cancellation depends in turn upon a clean separation of the time scales, 

1. r qu , of the very rapidly oscillating modes of the microscopic quantum theory, 

2. t c i, of the more slowly varying mean number of particles in the adiabatic number basis, 
and 

3. t p i, of the collective plasma oscillations of the electric current and mean electric field pro- 
duced by those particles. 

In the limit r qu <C r c i quantum coherence between the created pairs can be neglected because of 
efficient dephasing and a (semi) classical local kinetic approximation to the underlying quantum 
theory becomes possible. In the limit t c i <C t p i the electric field may be treated as approximately 
constant over the interval of particle creation. Thus when both inequalities apply we can replace 
the true nonlocal source term which describes particle creation in field theory by one that depends 
only on the instantaneous value of the quasistationary electric field, at least over very long 
intervals of time. 

The essential mathematical ingredient in the exploitation of this hierarchy of time scales is 
an asymptotic expansion of the wave functions and particle number for constant electric fields 
uniformly valid on the real time axis, so that secular particle creation effects (which are lost in the 
usual nonuniform WKB expansion) are retained. It is this precise sense of evaluating the effect of 
rapid degrees of freedom on slow degrees of freedom by treating the latter as constant in leading 
order of a uniform asymptotic expansion (which recalls the Born-Oppenheimer approximation 
in atomic and molecular physics) and deriving a local effective source term for the change of 
adiabatic particle number that we refer to as the Markov limit of the quantum Vlasov equation. 

The importance of a uniform asymptotic expansion of the wave functions is that secular 
particle creation effects are retained in an expansion valid everywhere on the real time axis. The 
true wave functions exhibit a sharp change in amplitude, on the time scale t c i, at or near the time 
of the semiclassical creation event which is captured very well by a uniform asymptotic expansion 
in terms of Airy functions. As we shall see, if one is interested only in the collective phenemena 
on time scales of t p i or longer, then the details of the particle creation process on the time scale 
t c i are unimportant and one can replace the momentum distribution of the source term by one 
localized at zero kinetic momentum, as has been the practice in the earlier phenomenological 
approaches, provided only that the integrated distribution gives the correct total creation rate. 
This will clarify the precise conditions of validity of such anatze for the first time. 
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Since an asymptotic (not a convergent) expansion is involved, the limit of the ratio of time 
scales r qu /r c i —>■ for fixed t and the long time limit t — > oo of the evolution for fixed ratio r qu /r c i 
do not commute in general. Hence for any small but finite ratio r qu /r c i there can be eventually a 
very large but finite t at which the quantum phases reassemble and the irreversible local kinetic 
description breaks down. Up to this very long (typically exponentially and possibly infinite) 
recurrence time the system behaves in many practical respects like an irreversible one, in which 
the quantum phase coherence between the created pairs appears to have been lost. In this way 
the apparent incongruity of an effectively irreversible time evolution emerging from a unitary 
Hamiltonian field theory is removed. 

The paper is organized as follows. In the next section we review (scalar) QED mean field 
theory in the leading order of the large N expansion. By exhibiting explicitly the Hamiltonian 
structure of these equations we demonstrate that they are completely time reversible. In Section 
3 we define the adiabatic particle number basis which is selected by the Hamiltonian evolution 
and derive the exact nonlocal form of the quantum Vlasov equation for this quantity. The 
quantum density matrix in this basis is also derived. In Section 4 we solve for the source term 
of the Vlasov equation in the limit of constant mean electric field, and study the pair creation 
process for this case in some detail. It is shown that particle creation in a fixed external field 
produces an electric current which grows linearly with time, so that any amount of particle 
creation (no matter how small) eventually requires a substantial backreaction on the field in any 
self-consistent treatment. In Section 5 the technique of uniform asymptotic expansions for the 
mode functions and adiabatic particle number is brought to bear. The source term for particle 
creation in a constant field is calculated to leading order in this asymptotic expansion in terms 
of Airy functions and yields an effectively Markovian source term for the local Vlasov equation 
describing pair creation in weak, slowly varying electric fields. The circumstances under which 
further approximation of the Airy function source term by an instantaneous delta function source 
term becomes permissable is also discussed. In Section 6 the dynamical backreaction problem 
for the charged particles whose current is self-consistently coupled to the mean electric field is 
compared to the two (Airy and delta function) local approximations for the Vlasov source term 
in the kinetic description, and relatively good agreement is obtained. We close with a summary 
of our results and some concluding remarks on possible generalizations of the analysis to other 
systems of interest. The derivation of the density matrix in the adiabatic particle number basis 
is relegated to an Appendix. 



II. SCALAR QED IN THE LARGE N LIMIT 

Let us begin by reviewing the equations of motion for scalar QED in a uniform electric field 
in the semiclassical limit in which the matter field is fully quantized and the electromagnetic 
field is treated classically. This limit can be obtained in a consistent way by taking the leading 
order of a large N expansion (where N is the number of identical copies of the charged matter 
field) We take the electric field spatially homogeneous, and express the vector potential in 

the gauge, 

A = A(t)z, A = 0, (2.1) 

so that the electric field is 

E = -Az = Ez. (2.2) 
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The charged scalar field operator is expanded in Fourier modes in Fock space in the usual way, 
$(x, t) = ±= £ e* k > k (t) = E {e lk ' x /k(t)a k + e-^f\{t)b{) . (2.3) 



V k k 
The time-independent creation and destruction operators obey the commutation relations 

[flk, 4'] = [ fe k, 4] = ^k,k' (2.4) 

in the finite large volume V, and the Fourier components, 

Vk(t) = /k(t)a k + /i(t)6l k (2.5) 

may be regarded as (complex) generalized coordinates of the field $ for the purposes of the 
Hamiltonian description. The momentum canonically conjugate to this coordinate is 

7r k (t) = 4(t) = fc(t)at + / k (t)6_ k , (2.6) 

which obeys the canonical commutation relation, 

[<Pk, Tk'] = «^k,k' , (2.7) 

provided that the mode functions satisfy the Wronskian condition, 

/k/ k - A/ k = ih , (2.8) 



and (|2.4j) is used. 

The time dependence in this basis is carried by the complex mode functions / k (i) which 
satisfy the equations of motion, 

+<«£(*) W) = (2-9) 

where the time dependent frequency cu k (t) is given by 

ujl(t) = (k - eA) 2 + m 2 = (k - eA{t)f + k 2 ± + m 2 . (2.10) 

Here is the constant canonical momentum in the z direction which should be clearly distin- 
guished from the gauge-invariant but time dependent kinetic momentum, 

p(t) = k-eA{t) ; p = -eA = eE (2.11) 

which reflects the acceleration of the charged particle due to the electric field. In the directions 
transverse to the electric field the kinetic and canonical momenta are the same and do not need 
to be distinguished, i.e. we shall use the notation p± = k± interchangeably. When expressed as 



a function of the kinetic momenta we use the notation, u(p, p±) = \jp 2 +p]_+ rn 2 , or simply uj. 
The mean value of electromagnetic current in the z direction is 

j(t) =2eJ [dk] (k - eA{t)) |/ k (t)| 2 (l + iV + (k) + AL(-k)) (2.12) 



where 



N + (k) = (4a k > 

AL(k) = (blb k ) (2.13) 
are the mean numbers of particles and antiparticles in the time independent basis and 

^j^Swf (2 ' 14) 

in the infinite volume continuum limit. We make use of the freedom in defining the initial phases 
of the mode functions to set the correlation densities (dkCik) = (&k&k) = 0, without any loss of 
generality. 

The mean charge density must vanish, 

f(t) = e J [dk] [N + (k) - AL(-k)] = , (2.15) 

by Gauss' Law for a spatially homogeneous electric field (ie.V • E = 0). We shall further restrict 
ourselves to the subspace of states for which 

iV+(k) = AL(-k) = iV k (2.16) 

for simplicity in what follows, although this is a stronger condition than is required by Eq. (|2.15|) . 
Clearly the vacuum iV+(k) = iV_(— k) = (as well as a thermal mixed state) belongs to this 
class of states. 

Self-consistent evolution of the mean electric field requires coupling it to the expectation 
value of the current of the charged field by the only nontrivial Maxwell equation remaining in 
this homogeneous example, namely, 

- E = A = j = 2eJ [dk] (k - eA(t)) |/ k (t)| 2 (l + 2iV k ) . (2.17) 

For the analysis of the source term in a constant electric field and its uniform expansion in the 
next three sections we will treat the electric field as fixed and nondynamical, returning to Eqn. 
(12.171 ) and the dynamical backreaction problem in Section VI. 



By a slight change of notation it is possible to recast the mean field evolution equations fl2.9| ) 
and (|2.17|) together with the quantum Wronskian condition ( |2.8|) as Hamilton's equation for an 
effective classical Hamiltonian in which h appears as a parameter. Defining the real quantities, 

a k = 1 + JV+(k) + AL(-k) = 1 + 2iV k , 
= ^k|/k(t)| 2 , and 
7fr(t) = Ut) (2.18) 

we find that the mode equation Q2.9|) can be rewritten in the form, 

4 = ^k = -^k + ^, (2.19) 



when account is taken of (|2.8| ). This last equation together with the Maxwell equation ( 2.17|) 
will be recognized as Hamilton's equations for the Hamiltonian, 
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H eff (A,p A ; {£ k }, fo}; {a k }) = ^ + E (vl + ofel + ^) , (2-20) 

where pa = —E is the momentum conjugate to A and 7/ k is the momentum conjugate to £ k - 

Moreover, the quantum statistical density matrix of the charged scalar field corresponding to 
the mean field evolution can be written as a product of Gaussians in Fourier space, viz. 

(KIIpIW) = n<K}|p(£k,^ k )|{^k}> = Up^ (2-2i) 

k k 

with 

Pk = (2< k ) ^ exp j [<p k <p k + <p k <p k ] + ? — [^ k ^ k - <pk<pk] + [<p k <p k + ^ k ]| , 



(2.22) 



and <p k is the complex generalized coordinate of the classical field in Fourier space, defined 
by Q2.5p (with a k and 6_ k treated as c-numbers). The Liouville equation for the evolution of 
this density matrix according to the quantum Hamiltonian of a free charged scalar field in a 
background electric potential, 

P = -i[Hgu, p}; H qu = - (^4 + uil^vi + h.c.) (2.23) 

1 k 

gives precisely the equations of motion ( |2.19| ) for the width parameters of the time dependent 
Gaussian. The effective classical Hamiltonian ( [2.20| ) is nothing else than the expectation value 
of the quantum Hamiltonian of scalar QED H qu in the Gaussian density matrix p, i.e. H e ff = 
Tr(pH qu ). Notice that in this Schrodinger representation of the time evolution all the equations 
are local in time, i.e., they involve a single time argument, and there is no need to introduce 
Wigner functions with two time arguments, although these correlation functions at unequal times 
may be calculated easily enough from knowledge of the density matrix, if desired. In contrast 



to several earlier approaches to kinetic theory from field theory principles ||r9| , |20| , we shall not 
require these unequal time correlators or Wigner functions. 

The constant parameters cr k = 1 + 2iV k > 1 measure the extent that the quantum state is a 
mixed state. If iV k = 0, a k = 1 and the state is pure, as is evident from the vanishing of the last 



term in (p.22j), so that the density matrix becomes a simple product |^)(^|. In either the pure 



or more general mixed state the density matrix ( 2.22| ) possesses a U(l) symmetry under 



(pt -> y? k exp« k ) 

<Pk ex P«k) (2.24) 

for each k. This is a reflection of the fact that the generator of the local U(l) gauge transformation 
of electrodynamics for a spatially uniform electric field is the charge density ( 2.15| ), and we have 



restricted ourselves to charge symmetric states obeying fl2.16Q , so that the density matrix has 
this U(l) invariance in each Fourier mode independently. 

In this leading order of the large N expansion the density matrix of the electric field is also 
a Gaussian and multiplies the matter field Gaussian above, so that the evolution of the closed 
system with the backreaction Eqn. (|2.17 ) is also Hamiltonian. Clearly the Hamiltonian evolution 
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equations ( |2.9| ), with or without the Maxwell Eqn. ( |2.17| ) are completely time reversible upon 
reversing the signs of all the momenta. 

Forgetting for the moment the Maxwell equation of backreaction on the electric field we see 
that the mean field evolution is equivalent to a set of time- dependent harmonic oscillators, with 
a different time-dependent frequency w^t) for each Fourier mode f^. Treating these frequencies 
as arbitrary, slowly varying functions of time we may write down the Hamilton- Jacobi equation 
corresponding to the effective classical Hamiltonian H e ff, namely, 



'dW k 



,2 2 



(2.25) 



and find that the Hamilton principal function evaluated over one full period, 



1 



2nh 2irh 



\ 



2e k - u££ - 



h 2 aj 
4£ 



Ok 

2 



(2.26) 



is an adiabatic invariant of the periodic motion. Since <7k is strictly a constant for all k, this 
implies that 



2JV k (i) 



(2.27) 



is an adiabatic invariant of the motion for slowly varying ouk(t). It is this adiabatic invariant that 
defines a time-dependent particle number basis which becomes the appropriate one for making 
contact with the Boltzmann-Vlasov kinetic description of particle creation. 



III. THE ADIABATIC NUMBER BASIS 



From the field theory development of the last section we note that N +} N_ and k, appear 
quite naturally in either the time independent (Heisenberg) or time dependent (Schrodinger) 
descriptions as constants of motion under the Hamiltonian evolution. However, kinetic theory 
is expressed in terms of time-evolving quantities J\f + (t),J\f_(t) and p(t) which must be clearly 
distinguished from the analogous time independent quantities above. The difference between the 
canonical and kinetic momenta k and p(t) in (|2.11|) is clear enough on basic kinematic grounds. 
The specification of the time dependent particle numbers Af+(t) and A/l(t) may not be quite 
as obvious, but as they provide the essential connection between the field theory and kinetic 
descriptions we must take special care to be equally clear and explicit about their definition. 
This requires that we introduce a Bogoliubov transformation from the time independent to a 
time-dependent (but adiabatic) number basis. 

The observation underlying the introduction of this basis is that the mode equation ( [2.9|) 
generally posseses time dependent solutions which have no clear a priori physical meaning in 
terms of particles or antiparticles. The familiar notion that positive energy solutions to the wave 
equation correspond to particles while negative energy solutions correspond to antiparticles is 
quite meaningless in time dependent background fields where the energy of individual parti- 
cle/antiparticle modes is not conserved, and no such neat invariant separation into positive and 
negative energy solutions of the wave equation is possible. This is just a reflection of the fact 
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that physical particle number does not correspond to a sharp operator which commutes with 
the Hamiltonian, i.e. particle/antiparticle pairs are created or destroyed, and physical particle 
number is not conserved in time dependent background fields. 

Given this fact, one possible point of view is to forget completely about particle number in 
time-dependent backgrounds and deal only with conserved physical currents like j(t) in ( J2.12 ). 



Indeed, in arbitrarily strong and rapidly time-varying fields this is the only possible point of 
view, since all notion of even an approximately conserved particle number disappears, and there 
is no possibility whatsoever of a classical kinetic description in such extreme situations. One 
must rely then exclusively on the field theoretic framework. 

When the fields are not quite so strong and/or so rapidly varying in time we would expect to 
be able to define a particle number which varies slowly enough for the comparison to an effective 
semiclassical kinetic description to be meaningful. Clearly this physical slowly varying particle 
number is not the iV k of the time independent Heisenberg basis defined by ( [2.13| ) above, since this 
iV k is part of the initial data, a strict constant of the equations of motion, no matter how strong or 
rapidly varying the electric field is. The physical particle number at time t must be defined instead 
with respect to a time dependent basis that provides some criterion to distinguish particles from 
antiparticles at the instant t. This time dependent basis (the adiabatic number basis) permits a 
semiclassical correspondence limit to ordinary positive energy plane wave solutions in the limit 
of slowly varying u>k(t), and is related to the Heisenberg basis by a time dependent Bogoliubov 
transformation. Since the mode equation (|2.9|) is the equation of motion of a (complex) harmonic 
oscillator with time varying frequency cj k (t), governed by the effective classical Hamiltonian H e ff 
of (|2.20| ) standard arguments from classical Hamilton- Jacobi theory inform us that there is an 



adiabatic invariant proportional to the energy of the oscillator divided by its frequency, and given 
by ( |2.27|) . It is this quantity which can be used to define an adiabatic particle number and to 
make the connection with classical kinetic theory. Corresponding to this slowly varying action 
variable there is a conjugate angle variable which is rapidly varying, of which classical kinetic 
theory takes no account. 

The adiabatic basis is defined by first constructing the adiabatic mode functions, 



/ k (t) = J — ^- exp (-% f ^(t')dt' ) . (3.1) 



2wk(t) 

We will make use of the shorthand notation for the phase, 

6 k (t) ee J* \o k (t')dt' , (3.2) 

suppressing the explicit dependence on t (and occasionally also the momentum index k) except 
when needed for clarity in most of the following. The lower limit of the integral in (|3.2j) and 
therefore also the absolute phase of the mode function / are left arbitary for the moment, to 
be fixed in a convenient way in the next section. It is clear that in the limit of arbitrarily weak 
electric fields u; k (£) becomes independent of time and can be removed from the integral in Q3.2| ). 
In that limit the adiabatic mode function becomes the usual positive energy plane wave solution 
with respect to which the usual definition of particle number is taken. Otherwise the adiabatic 
mode functions ( |3.1| ) will not be exact solutions of the mode equation ( |2.9| ), but we are still free 
to specify a basis with respect to them, provided only that w k (t) remains real and positive for 
all k and t. 
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The transformation to this basis from the original one is specified by the two linear relations, 

/k(t) = ak(t)/k(t)+/8k(t)/i(t) 

f k (t) = -iu k a k {t)f k (t) + iu k p k (t)f£(t) (3.3) 

between the exact and adiabatic mode functions. When the phase of / is fixed these relations 
completely fix the complex coefficients a k (t) and f3 k (t). It is straightforward to solve for the 
Bogoliubov coefficients directly in the form, 

«k = i(/k - «*Jk/k)/k and 

Ac = -«(/k + «^k/k)/k. (3.4) 

An equivalent form of this Bogoliubov transformation in the Fock space of creation and destruc- 
tion operators is 

a k = a* k {t)a k {t) - f3* k {t)b ] _ k {t) 

6l k = a k (t)bl k (t) - p k (t)a k (t) , (3.5) 

so that the field coordinate (p k (t) may be expressed equally well in the time independent basis 
by fl2.5|) or in the time dependent basis by 

Mt) = /k(t)a k (t) + /i(t)6l k (t) , (3.6) 

and likewise the field momentum variable is given either by ( |2.6| ) or by 

7r k (t) = -iu k (t)f£(t)a{{t) + iu k {t)f k {t)b„ k {t) . (3.7) 

The transformation from the time independent (a k , 6L k ) basis to the time dependent adiabatic 
basis (a k , 6L k ) requires two independent relations (|3.3|) or (|3.5| ), corresponding to a canonical 
transformation in a two dimensional (complex) phase space, for which 

KI 2 - IA<r = i (3-8) 

for each k. It is easily verified that (|3.4|) satisfies this relation when the Wronskian condition 
( |2.8| ) is used. Because of ( |3~8| ) the magnitude of the Bogoliubov transformation to the adiabatic 
number basis 7 k (t) may be specified by 

\a k (t)\ = cosh7 k (t) , 

\/3 k (t)\ =sinh 7k (t). (3.9) 
We now define the adiabatic particle number to be 

A/" k (t) = (al(t)a k (t)) 

= |a k | 2 (4 a k) + |Ac| 2 (fr- k fr- k ) 

= (l + |/? k | 2 )iV + (k) + |/? k | 2 (l + iV_(-k)) 

= N k + (1 + 2iV k ) \(3 k (t)\ 2 

= N k + (1 + 2N k ) sinh 2 7k (t) . (3.10) 
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The second of the relations ( |3.3| ) is essential to define the adiabatic basis in which particle number 
is given by the ratio of energy to frequency. In fact, 



hoj-kif) hook 
= (l + 2JV k ) (l + 2|/5 k | 2X 
= l + 2jV k (i). (3.11) 

Hence the particle number W k (/!;), though time dependent, is an adiabatic invariant of the 
motion. Consequently, it is the natural candidate for a particle density in phase space for a 
kinetic description, becoming the ordinary asymptotic constant particle number in the limit of 
slowly varying u; k (t). Although this choice of basis is not unique, since we could have chosen a 
different condition on / k in ( |3.3| ), it is the only basis where the ratio £ k /u; k is simply related to 
particle number (without the appearance of a> k or higher derivative terms, for example) which is 
the standard adiabatic invariant of the harmonic oscillator with time dependent frequency as in 



( 2.27 ). In different contexts (such as particle creation in external gravitational fields) it may be 
appropriate to consider slightly different definitions of the adiabatic number basis. 

Now that we have completely specified the time dependent particle number basis it is straight- 
forward to derive the equation of motion which it obeys. We note that from the explicit repre- 
sentation ( |3.4| ) by differentiation and use of the mode equation ( |2.9| ) we have 

d k = ^/? k exp(2i6 k ) , and 

4 = 7^ L akexp(-2z6 k ). (3.12) 
2^ k 

These two first order differential equations are entirely equivalent to the second order mode 
equation in Hamiltonian form. We now obtain by differentiating (|3.10|) , 

^-AT k = 2(l + 2iV k )Re(^ k ) 
at 

= ^ (i + 2N k ) Re{a k /? k exp(-2*0 k )} 

= — Re{C k exp(-2i6 k )} , (3.13) 

where we have defined the time dependent pair correlation function, 

C k (t) = (a k (t)6_ k (t)) 

= (l + 2iV k )a k /3*. (3.14) 

Thus the time derivative of the adiabatically slowly varying particle number involves the pair 
correlation function C k (t) which is itself very rapidly varying, since the time dependent phases 
on the right side of (|3.14j ) add rather than cancel, although the phases do nearly cancel in the 



final combination of ( |3.13| ). The time derivative of the pair correlation function, 



d C k = p~ (1 + 2iV k ) exp(2*6 k ) (l + 2|/? k | 



dt 2u; k 
d; k 



} ._ + 2AT k )exp(2^e k ), (3.15) 
2^ k 
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brings us back again to Ak- This last equation may be solved formally for C k and subsituted 
into (|3.13|) to obtain 

4;A/k = ^-f dt' {^(0 (1 + 2jV k (t')) cos [26 k (t) - 26 k (0]| , (3.16) 

where we have assumed that C k vanishes at some t = to (which could be taken to — oo). 

Equation fl3.16| ) may be be called a "quantum Vlasov equation," in the sense that it gives 



the quantum creation rate of particle number in an arbitary time varying mean field. Let us 
remark that the Bose enhancement factor (1 + 2jV k ) appears in ( 3.16Q , so that both spontaneous 



and induced particle creation are included automatically in the quantum treatment. The most 
important feature of Eqn. ( |3.16| ) for our present purpose is that it is nonlocal in time, the particle 



creation rate depending on the entire previous history of the system. In that sense the particle 
creation process is certainly non-Markovian in general [T5|Jl7| . Eqn. ( |3.16|) becomes exact in 



the limit in which the electric field can be treated classically, i.e. the large N limit in which 
real and virtual photon emission is neglected, and there is no scattering. Inclusion of scattering 
processes lead to collision terms on the right side of ( p.!6|) which are also nonlocal in general. This 
nonlocality is essential to the quantum description in which phase information is retained for all 
times. The phase oscillations in the cosine term are a result of the quantum coherence between 
the created pairs, which must be present in principle in any unitary evolution. However, precisely 
because these phase oscillations are so rapid it is clear that the integral in receives most of 

its contribution from t' close to t, which suggests that some local approximation to the integral 
should be possible, provided that we are not interested in resolving the short time structure or 
measuring the phase coherence effects. The time scale for these quantum phase coherence effects 
to wash out is the time scale of several oscillations of the phase factor @ k (£) — Qk(t'), which is 
of order r qu = 27r/c<j k = 27rfr/e k , where e k is the single particle energy. 

The steps we have just performed to arrive at (|3.16|) are a special case of the general projection 
formalism of Zwanzig |2T|, where some subset of fast dynamical variables deemed "irrelevant" 



(in this case Ck) are eliminated in favor of slow variables deemed "relevant" (in this case A/" k ) 
Because the two variables are coupled by the underlying Hamiltonian equations of motion the 
result of solving for some variables in terms of others is generally nonlocal in time. The nonlocal 
form ( |3.16| ) is still completely equivalent to the mode equation ( |2.9| ) and absolutely nothing has 
been lost (or gained) by this rewriting. In other words, the projection method is essentially 
free of any physical content, until and unless one makes further approximations that replace the 
nonlocal relations satisfied by the relevant observables by local ones. It is at this point that great 
care must be exercised, since the precise form of the local approximation made will determine 
the usefulness and range of validity of the resulting truncation. 

For example, one's first idea might be to remove the Bose enhancement factor l + 2A/'k(£ / ) from 
the integral, or simply ignore it entirely, on the basis that it is slowly varying function (for real 
t') , and attempt to perform the remaining integral over the rapidly varying phase by the method 
of stationary phase. However, this does not lead to a useful result, since the phase becomes 
stationary at G k = ^ k = which is precisely where the integrand has a pole in the complex 
t' plane. This is more than just an inconvenience since the vanishing of is the condition 
of a turning point, where the adiabatic approximation certainly breaks down. Hence it is in 
just the neighborhood of this point in the complex t' plane where the dominant contribution to 
the integral arises that the removal of the factor 1 + 2A/ k (£ / ) from the integrand or its neglect 
cannot be justified. The stationary phase method has been carried out nonetheless for certain 



12 



background field problems where integrals similar to that in ( 3.16 ) arise, with the result that the 
correct exponential factor but the incorrect prefactor is obtained ||. 

The importance of the complex turning point (s) for determing the asymptotic mixing between 
particle and antiparticle modes as t — > ±00 has been emphasized by Marinov and Popov in refs. 
0. In their method the analytic continuation of the solutions of the mode equation around 
the Stokes' lines emanating from the turning point in the complex time plane determines the 
sub dominant component of the wave function with the opposite sign of the frequency on the real 
axis. The amplitude of this exponentially subdominant component of antiparticle waves in the 
wavefunction is the Schwinger particle creation effect. However, the method outlined by these 
authors does not seem to be applicable to the integral in ( |3.16| ) directly, since it is designed 
for calculating the particle creation asymptotically over infinite time, not for determining the 
evolution of the particle creation process in finite real time t, which is what we require for the 
transport description. 

Finally, if one takes no account of the stationary phase point in the complex t' plane but 
attempts to approximate the integral in ( |3.16 ) entirely in real time, for example by integrating 
the rapidly varying cosine function by parts any number of times, it is easy to see that an 
asymptotic series is generated in which the exponentially small subdominant solution can never 
appear after any finite number of such steps. Any asymptotic expansion of the wavefunction 
on the real axis which discards the exponentially small antiparticle component will miss the 
Schwinger creation effect at late times. 

From this discussion we see that the essential difficulty with Eqn. ( |3.16| ) is that the point (s) 
in the complex t! plane where the phase 0k is stationary must play the critical role in determining 
the particle creation for asymptotically late times, but we cannot evaluate the contribution to 
the integral of these stationary phase points where vanishes without in effect knowing the 
full Ak, k>k an d @k as analytic functions in the entire complex t plane before we even begin. If 
we were in possession of these analytic functions we would already have the full solution to our 
dynamical problem, without any need to make any approximations to the integral. This is clearly 
impossible except for a small number of special cases where the complete analytic structure is 
known a priori. Thus the nonlocal form of the quantum Vlasov equation ( |3.16|) makes it difficult 
to extract any useful information about a source term for a kinetic description in general. 

Consideration of this difficulty immediately suggests a different approach. Instead of trying 
to work with the nonlocal equation ( |3.16| ), in the next section we evaluate the spontaneous pair 
creation rate ;§Ak(t) for a constant electric field analytically and directly in real time, thereby 
assuring agreement with the Schwinger result in both its exponential and nonexponential factors. 
This is one of special cases where A/"k and its time derivative can be evaluated analytically in local 
form, directly from the definition ( |3.10|) without any need for the nonlocal integral representation 
( 3.16|) . Then by making use of an asymptotic expansion of the exact analytic result for constant 
fields, uniformly valid everywhere on the real time axis, we obtain a useful local approximation 
to the spontaneous pair creation rate for the slowly varying electric fields, without any need 
for analytic continuation or stationary phase methods in complex time. By such an approach 
we shall bypass completely the difficulties of dealing with the nonlocal integral equation ( |3.16| ) 
resulting from the projection method. 

The transformation to the adiabatic number basis and elimination of the rapid variables 
Ck in favor of the slow variables Ak by Eqns. ( |3.14| )- (|3.16|) has its counterpart in the density 
matrix description as well. It is shown in the Appendix that the density matrix (|2.22|) may be 
transformed to the adiabatic number basis, with the general form of the nonvanishing matrix 
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elements given by ( |A27| ). In the pure state case <7k = 1 the only nonvanishing matrix elements of 
p are in uncharged pair states with equal numbers of positive and negative charges, £^ = = 
n k , with £k the number of pairs in the mode k, viz. 

(24|p|24) = e^'-^W sech 2 7k (t) (tanh 7k (t))' k ' +£k (3.17) 

£T=1 

where the magnitude of the Bogoliubov transformation, 7k(£) is defined by ( |3.9| ) and its phase, 
i?k(^) is specified by 

a k /? k e" 2i0k = -sinh 7k cosh 7k e i,?k . (3.18) 

Hence the off-diagonal matrix elements £' ^ £ of p are rapidly varying on the time scale r qu of 
the quantum mode functions, while the diagonal matrix elements £' = £ depend only on the 
adiabatic invariant average particle number via 



(24M24) 



= p 2£k = sech 7 k tanh k 7 k 



;i + i/5 k iT k+1 (i+A4r k+i 



(3.19) 



and are therefore much more slowly varying functions of time. The average number of positively 
charged particles (or negatively charged antiparticles) in this basis is given of course by 

oo 

E = A/k . (3.20) 

£ k =0 

Thus the diagonal and off-diagonal elements of the density matrix in the adiabatic particle num- 
ber basis stand in precisely the same relationship to each other and contain the same information 
as the particle number A/k and pair correlation Ck respectively. 

Using the representation ( 3.17|) or ( 3.19|) we can understand how entropy can increase and 



the evolution become time irreversible if we replace the exact nonlocal quantum Vlasov equation 
( |3.16| ) by a local expression in which the rapid phase variables Ck, $ k) or the off-diagonal matrix 



elements of p no longer appear. Time reversal in the field theory requires that both the slow and 
fast variables be time reversed, which involves the full density matrix p. If we restrict attention 
to only the diagonal matrix elements of p in the adiabatic particle number basis without any 
account of the phase information present in the rapidly varying off-diagonal elements, then time 
reversal no longer holds. In the effective density matrix ( |3.19| ) the diagonal elements p2e k may be 
interpreted (for <7k = 1) as the independent probabilities of creating 4 pairs of charged particles 
with canonical momentum k from the vacuum. This corresponds to disregarding the intricate 
quantum phase correlations between the created pairs in the unitary Hamiltonian evolution, and 
treating the creation events as essentially independent in a stochastic Markovian processes. Thus 
the Markov approximation to the field theory arises quite naturally when the quantum density 
matrix is expressed in the adiabatic particle number basis. 

Such an approximation is known to be quite accurate for long intervals of time in the back- 
reaction of the current on the electric field producing the pairs, for the simple reason that the 
phase information in the pair correlations cancels very efficiently when one considers the sum 
over all the k modes in the current (|2.12|) . It is for this reason that for practical purposes one can 
approximate the full Gaussian density matrix over large time intervals by its diagonal elements 
only, in this basis. Naturally this truncation of the unitary Hamiltonian evolution according to 
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(|2.23 ) leads to a nonunitary irreversible evolution in which the effective von Neumann entropy 
of the diagonal density matrix (|3.19|) , 



S eff{t) = -Trp e ff lnp e// = - Y, PW* ln ^ k (3.21) 

k £ k =0 

can increase with time. In fact, upon subsituting (|3.19|) , the sums over are geometric series 
which are easily performed, with the result that the von Neumann entropy of this truncated 
density matrix, 

Seff(t) = £{(l+A/k) ln(l + A/k) - A/k In A/k} (3.22) 
° =1 k 

is precisely equal to the Boltzmann entropy of the single particle distribution function A/k(t). 
Hence 



dt en Y V A/k /eft 

increases if the mean particle number increases. This is always the case on average if one starts 
with vacuum initial conditions, Ok = 1, since |/?k| 2 is necessarily nonnegative and can only 
increase if it is zero initially ||22|| . Locally, or once particles are present in the initial state, there 
is no reason why particle number or the entropy ( |3.23| ) must continue to increase monotonically 
in time, and indeed small temporary decreases are observed in backreaction simulations ||. 
Hence there is no Boltzmann H-theorem for the effective entropy ( |3.23| ) without introducing 
some explicit time averaging and/or further assumptions into the scheme. 

Before closing this section we wish to take note of one additional especially simple property 
of the adiabatic particle number basis. Inserting the Bogoliubov transformation of the mode 
functions ( |3.3|) into the expression for the current ( | . 1 2 ) we obtain 

j(t) = e [ [dk) {k ~ e ^ {t)) (1 + 2|/3 k (t)| 2 + 2Re{a k /?ke- 2j0kW })(l + 2iV k ) . (3.24) 

We note that the vacuum term in this expression, / [rfk] ^^(t) vanishes by charge conjugation 
symmetry, when proper gauge invariant integration boundaries are chosen. Using the mean value 
of particles in the adiabatic number basis ( 3.1 0|) , its time derivative and the equations of motion 
(|3.13| ), we can rewrite the current fl3.24|) as 

j(t) =2eJ Mk] (fc ~^ (t)) A/k(t) + | / [dk]^ k (t)A/ k (t) 

= jcond "I" jpol ■ (3.25) 

On the other hand, from a classical point of view if the particle distribution A/k is coupled to a 
uniform electric field the energy density and its time derivative are given by 

E 2 f 

e = — + 2 J [dk]cjkA/"k (3.26) 
i = EE + +2 J [dk] ^eE ^~ eA V k + uj^^j = . (3.27) 
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Using the Maxwell equation —E=j this last relation is precisely the same as the mean value 
of the quantum current in ( |3.25| ). Hence we may identify the adiabatic particle number Ak(t) 
with the (quasi) classical single particle distribution. Other definitions of time varying particle 
number, such as that used in our own earlier work |l|] do not have this property or admit this 
simple quasiclassical interpretation. This exercise also demonstrates that the two terms in the 
mean current Q3.25|) should indeed be interpreted as the conduction and polarization terms of 
the earlier phenomenological descriptions. 



IV. CONSTANT ELECTRIC FIELD 

In order to derive the source term due to particle creation in a slowly varying electric field, 
we first analyze the time structure of the creation process in a constant, uniform electric field, 
for which 

A(t) = -Et. (4.1) 

It is useful to define the rescaled dimensionless variables 

k + eEt ep(t) B + m 2 

u = e . = ~p=L and A = \ - > , 4.2 

^El ,^E\ \eE\ 1 } 



where e = e(eE) = ±1 is the sign of eE. Then the mode equation ( |2.9|) may be put into the 
form, 

(^ + " 2 + A ) / = <«> 

whose solutions are parabolic cylinder (Weber) functions. In fact, the two complex conjugate 
pairs of solutions, 

/(+)(«) = /(-)(«) « D -\+i\ H 1 ~ and 
f(+)(u) = oc D_i_ ih ((1 + i)u) (4.4) 

each comprise complete sets of basis functions in which to expand the scalar charged field $. 
Normalizing these solutions according to the Wronskian condition (2~B) and defining the phase, 



^ = - - -lnA + -ln2- - (4.5) 
^ 4 4 4 8 V ; 

we can write the properly normalized positive frequency mode functions in the form, 

/(+)k(t) = \leE\-^e~^e^ D_x +i x (-(1 - i)u) and 
fl+\t) = \2eE\-^e-^e-^ D_^_q ((1 + i)u) (4.6) 

which approach the adiabatic functions fk(t) in the asymptotic limits t — > — oo and t — >• oo 
respectively. Notice that with u defined including the e function as in ( |4.2| ) these limits are 
equivalent to u — >• — oo and u — >• oo respectively, independently of the sign of eE. The complex 
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conjugates of these solutions are the corresponding negative frequency mode functions and are 
denoted by /(-)k or / k ~^ respectively. The phase ip has been defined in such a way that the phase 
of the exact mode functions /k agrees with the adiabatic mode functions fl3.1| ) with phase 0k 
measured from the symmetric point u — 0, i.e. 



e k (t) 



u=Q 



dt'uk{t') 



du'Vu' 2 + A 



^uVu 2 + A + ^ In ^ 



it 



(4.7) 



It may seem surprising at first sight that the exact mode functions approach the adiabatic 
ones in the infinite past and infinite future even though the electric field E is constant and never 
vanishes in these limits. The reason for this is that the corrections to the lowest-order adiabatic 
mode functions involve 



5ul(t) 1 a) k 



0J\ 



2ui 



3u 2 



4^ k 



-2»u 2 + 2A) 
4{u 2 + A) 3 



(4- 



which goes to zero like |t|~ 4 as t — > ±oo. 

If the state of the system is the vacuum in state then the mode function /k to be used in 
is the /(+)k °f ( |4.6D and the effective source term for the creation of particles from the vacuum is 



7V k =0 



\%eE\ 



i 

"2e" 



4 



d_ 

dt 



' d 



D 



4 



l)U) 



(4.9) 



We note that for a strictly constant electric field this is an exact result for the rate of adiabatic 
particle number change starting from vacuum initial conditions at t — — oo. Phase correlation 
information for this particular initial state has not been discarded, although the pair correlation 
function does not appear explicitly in ( |4.9| ) which is local in time. 

Now since the two pairs of complex functions /(±)k and both satisfy the same second order 
wave equation there exist linear relations between them. Indeed it follows from the properties of 
the Weber parabolic cylinder functions that 



/ (+)k = afi +) + Pft ] 



(4.10) 



with mn 



a: 



Z7T 



r(i^) 



-ie 



2 



tt,\ 
4 



and 



(4.11) 



The fact that (3 ^ is the statement that the Bogoliubov transformation between the two 
basis pairs is nontrivial and the adiabatic vacuum state in the infinite past contains particle- 
antiparticle pairs with respect to the adiabatic vacuum state in the infinite future. The magnitude 
of this total Bogoliubov transformation from t = — oo to t = +oo is finite and given by 



sinh 2 7 



-VTA 



(4.12) 
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which is independent of k in the direction of the electric field. 

By transforming the Gaussian density matrix corresponding to the evolution of the charged 
scalar field in a background electric field one can show that ( |4.12| ) is also the mean number of 
particles in the final state with respect to the out vacuum, assuming that the field was prepared 
in the in vacuum. The details of this transformation are given in the Appendix. From the result, 
Q3.19Q with 7k replaced by 7 and the discussion of the previous section discarding the rapidly 



varying off-diagonal elements of p, we may interpret the diagonal elements as the probability of 
finding I pairs at late times if none were present initially. Hence the I = matrix element, 

sech 2 7 = (1 + e- nX y 1 (4.13) 

is the probability of creating no pairs in the given mode, and the probability that the vacuum 
remains the vacuum in the future is given by the product over all modes, 

J](l + e-^y 1 = exp (- J2 ln(l + e~" A )^ . (4.14) 

Taking the infinite volume limit this can be expressed as exp(— VTT) where the rate of vacuum 
decay per unit volume is 

r = 1 J[dk] ln(l + e- A ) . (4.15) 

Since the kinetic momentum of the created charged particles in the direction of the electric field 
is k + eEt, the longitudinal integration element dk can be replaced by eET as T — > 00 and the 
vacuum decay rate becomes 



eE 
eE 



J d 2 k ± ln(l + e- nX ) 

r 00 
/ d 2 k ± £ 



' g— nnX 



(2ir)»J n 

which is Schwinger's result for scalar QED. 

One should note that the replacement of the longitudinal momentum integral over k in (|4.15|) 



by eET in the large T limit can be justified only if one understands the time evolution of the pair 
creation event, for otherwise the expression ( |4.15 ) is formally meaningless. This replacement of 



/ dk by eET and the resulting finite expression (|4.16| ), which can be obtained by quite different 



methods imply that only those k in a linearly growing window in time actually contribute to the 
rate, although the mixing coefficient j3 over all time is independent of k. It is the time dependent 
evolution of ^(t) which we can investigate in detail with our definition of the time dependent 
adiabatic number basis in the next section. This definition smoothly interpolates between the 
in and out vacuum states specified respectively by the two wave functions in ( |4.6|) , so that /3k(t) 
starts at zero as t — > —00 and must approaches (3 as t — > +00. The wave functions depend on k 
and t only through the variable u defined in ( |4.2|) , and the potential u 2 + A is even in u. Hence 
we should expect each k mode to go through its creation event at a different time t according 
to k + eEt ~ 0, i.e. for the particles to be created with kinetic momenta near zero. We shall 
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see that this is indeed the case and that therefore the range of k which have gone through the 
creation process at time t depends linearly on t, which justifies the passage from ( [4.15 ) to ([4.16 ). 



Omitting the integration over kj_ and the phase space factor 1/(2-71") in ( |4.1(j| ), one obtains 
the probability per unit time per unit volume to produce pairs with transverse momentum 
]5,n]|. This result has been interpreted as the rate at which pairs are created 
a source term in Vlasov equation, which involves particle production |12 



m 



and used as 
However, a necessary 

condition for this interpretation to be correct is that the time integration over the rate of particle 
production 



eE 



dtln(l 



-ttX\ 



(4.17) 



be identical to the total number of particles produced per unit volume with transverse momentum 
k^, which is given by integration over k of Eq. ( f4.11| ) 



dk e 



-tt\ 



eE dte 



(4.18) 



Expressions ( [4.17]j and ( |4.18| ) are not equivalent, because the probability rate of particle produc- 
tion differs from the production rate of the mean value of particles. They become equal only in 
the limit of large A when both the probability and mean number of produced particles become 
very small. It is clear that the source term for the mean rate of particle production in the Vlasov 
equation should involve the latter quantity ( 4.18 ) in principle, without the appearance of any 
logarithm in the final answer. 



V. UNIFORM ASYMPTOTIC EXPANSION OF THE SOURCE TERM 

Eqn. ( [4.9| ) is the source term due to particle creation, with a specific choice of initial conditions 
and phase correlations in the initial state (namely none). Since all quantities in ([4.9|) are local 
functions of time, specified in terms of the mode functions ( |4.(j| ) there is no need to resort to 
the nonlocal integral equation (|3.16| ), and the time evolving phase correlation Ck need not be 
considered explicitly. Because of Eqns. (|4 . 1 0|) and ( |4. 1 1|) the Schwinger pair creation amplitude is 



certainly contained in (|4.9| ). Since our objective is the derivation of an effective Markovian source 
term for the Boltzmann- Vlasov equation fields which are slowly varying in time we now introduce 
the second important ingredient in our approach, i.e. the uniform asymptotic expansion of ( [4.9| ) 
for weak and slowly varying electric fields. 

In order to motivate the introduction of this asymptotic expansion observe that for a constant 
electric field each time derivative of the mode function brings with it a factor of 1/A. This 
can be made explicit by introducing a rescaled variable v which is independent of the strength 
of the electric field, u = v V\ and rewriting the wave equation ([4.3|) in the form, 



+ f 2 + l / = 0. (5.1) 



A 2 dv 2 



Next, when we allow the electric field to vary in time we can consider the standard adiabatic 
expansion for the mode function in the time-varying field fj]J 
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exp i J dt' Qk(t')\ 
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(5.2) 

for which (|3.1| ) is the lowest order term, corresponding to no derivatives in (|5.2|) and order A in 
the constant field case. The quasistationary or adiabatic approximation to the mode equation is 
obtained by treating the derivative terms as small compared to the leading order term, i.e. 

— < 1 and — < 1 . (5.3) 

uo 6 oo z 

In the case of a general electric field this implies that the field is both slowly varying and weak. 
For a constant electric field the adiabatic condition reduces to A 3> 1. 

Iterating the expansion to adiabatic order q, terms with q time derivatives of the adiabatic 
frequency uo^ in the general time varying field will appear together with terms with q powers of 
1/A in the constant field case. Thus there is a one-to-one correspondence between the asymptotic 
expansion of |/9k| 2 and J^|Ac| 2 m powers of 1/A in a constant electric field background to a local 



expansion of the current (|3.24j ) in higher time derivatives of the electric field in the general 
case. Now the transport approximation amounts to a truncation of this expansion at the lowest 
order required for a consistent backreaction dynamics. This is determined by the order of the 
backreaction equation, j = —E = A, which is second order in time. Thus we should expand the 
particle number \/3\ 2 only to second order i.e. 1/A 2 , in order to match the asymptotic expansion 
of the current to the order of the backreaction equation for a weak, slowly varying electric field, 
self-consistently determined by solving the Maxwell- Vlasov system. To retain higher orders than 
this in the current would not be correct mathematically, since these higher orders would also 
involve higher derivatives of E in the general time varying electric field, and such terms can never 
be calculated correctly by the constant E approximation of ( f4.9| ). At adiabatic order 2 the only 
effect of approximating the source term for a slowly varying electric field by (|4.9|) , evaluated in 
a constant field is the absence of the E term in the current ( |2.12j ) generated by the adiabatic 
expansion ( |5.2|) . This term is responsible for charge renormalization in the mean field theory |l],|2|. 
Hence for comparison between the mean field evolution and that of the Vlasov-Maxwell system 
one must specify the scale of the renormalized charge of the mean field theory by some other 
criterion, or it will differ in general from the classical charge appearing in the Vlasov equation 
by a finite renormalization. This precise correspondence we fix by a linear response analysis in 
Section VI. 

Even if we could calculate higher order terms (by calculating the source term in some other 
time varying background, for example) it would not be correct physically to include them since 
they would change the order of the Maxwell equation E = —j by making j a function of higher 
derivatives of E and thereby introduce unphysical high frequency runaway solutions, not present 
in the underlying microscopic quantum field theory, in a manner similar to the higher derivative 
Lorentz radiation reaction force. The important physical point here is that the order of the 
backreaction equation for time-varying electric fields determines the order of the asymptotic 
expansion we should use for the current, in the limit of weak, slowly varying electric fields, which 
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is the only limit in which such a replacement in the current is justified. The fact that the leading 
order asymptotic expansion of the constant field adiabatic particle number is already 1/ A 2 (as we 
shall see shortly) which is the highest order we need to go in the expansion, justifies the use of 
the constant field expression Q4.9|) , evaluated to this asymptotic order, for the local source term 
in the Markov limit of the quantum Vlasov equation. If the higher order terms in the expansion 
are numerically significant, then that is the signal that we must abandon the Boltzmann- Vlasov 
description entirely and return to the underlying field theory without the possibility of making 
any simple transport approximation to the self-consistent backreaction problem. 

The key point is that we require an asymptotic expansion of the mode functions and adiabatic 
particle number source term in ( |4.9|) in powers of 1/A that is uniformly valid in time (and 
longitudinal momentum k), in order that the exponentially small Schwinger amplitude (3 which 
is the only secular effect of particle creation which survives as t —>■ oo will not be lost in the 
expansion. This condition is not satisfied by the naive asymptotic expansion of / in simple 
exponential functions such as fl5.2| ). What apparently has not been so generally well recognized 
is that this failure of the usual adiabatic expansion to capture exponentially small (but secular) 
particle creation effects is due to the nonuniformity of the naive asymptotic expansion with 
respect to the limits t — > ±00. This limitation can be removed by an asymptotic approximation 
uniformly valid everywhere on the real time axis. 

In the case at hand, the asymptotic expansion of the solutions of Eqn. ( |4.3| ) uniformly valid 
everywhere on the real time axis have been given by Olver [23]. Converting to the notations of 
the present paper, Olver's result may be written in the form, 



/(+)k(t) 




8=1 



Q-2 



(5.4) 



where the coefficients P 2 s and Q2S+1 are certain functions of v = -7= given by 
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the complex variables £(i>), w and z are defined by 
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(5.5) 



and 



w = — A£ = —z* , 

where z is defined in the plane cut along the positive imaginary axis by 



(5.6) 
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(5.7) 



and the 7 S are the numerical constants, 
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With the definition of the phase of z according to ( p.7| ) the complex argument of the Airy 
functions in ( |5.4|) varies along the contour depicted in Fig. 1, as the real time t or u ranges from 
— oo to +00. 




FIG. 1. The contour in the complex z plane along which the argument of the Airy functions in 
(5.4), (|5.9| ), ( |5.11| ), and ( |5.12 ) are to evaluated. The cut of the | root appearing in fl5,7| ) is taken along 
the positive imaginary z axis from to ioo. Following (|5.6| ) the corresponding contour in the complex 
w plane is a straight horizontal line displaced from the real w axis into the lower half w plane by 
with Re w decreasing as t increases. 



The terms we have written here explicitly determine the uniform asymptotic expansion of 
/(+)k up to order 1/A 2 . Since we are interested only in the lowest nonvanishing order in the 
expansion we could retain only the lowest order term in ( |5.4| ), substitute it into ( fO|) to obtain the 
lowest order source term in the Vlasov equation directly. Some care is required in this procedure 
since the argument of the Airy functions depends on A through (|5.7|) and the equations of motion 



(|3.12p will not be satisfied unless both sides of the equation are expanded consistently to the 
same order in 1/A. For this reason it is useful to retain one higher order in the asymptotic 
expansion than would seem necessary at first sight, in order to have a nontrivial check on the 
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algebra via the equations of motion. 

The corresponding asymptotic expansion for the time derivative of the mode functions uni- 
formly valid on the real axis is: 
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where the coefficient functions are 
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(5.10) 



From these expressions the uniform asymptotic expansion for the time dependent Bogoliubov 
coefficient /3k(t) is easily computed from its definition in (|3.4|) , namely, 



/5 k ~ y^e 4 e 



-»e h 



1 + 5S 7 *(^ 



(5.11; 



Since we have shown by Eqn. (|3.11|) that the particle number | (3^ \ 2 is an adiabatic invariant to 
leading order in the time derivatives of the background, the lowest order A term in the asymptotic 
expansion must be absent from the particular linear combination in (|5.11|) . Indeed with s = 0, 
Vq = Qo = 1 and the symmetric linear combination of Airy functions, z^Ai(z) + z~^Ai'(z), is 
of order A -1 , as is verified explicitly in relations ( |5.20|) and ( |5.21[) below, by using ( |5.6|) and the 

2 

further asymptotic expansion of these functions for \z\ ~ A3 — > 00. Any other linear combination 
of the same functions, and in particular the antisymmetric combination, ziAi(z) — z~ 4 Ai'(z) is 
of order A . Anticipating this result and substituting (|5.5| ) and (|5.10|) into (|5.11| ), we obtain 
simply, 



z^Ai(z) + z 4 Ai'(2;) 



z*Ai{z) - z *Ai'(z) 



u 1 

r H 

(w 2 + A)2 3w 



(5.12) 



correct to the leading nonvanishing order, A -1 . Squaring ( |5.12|) and taking its time derivative 
gives the asymptotic approximation to the effective source term defined in (|4.9| ). Since w and z 
are functions of u and A (equivalently, v and A) which depend only on the kinetic momenta p(t) 
and p± through (|4.2| ), the effective source term (for vacuum intitial conditions at t — —00) may 
be written in the form, 



d 

S vac {p,P±;E) = eE—\(3(p,p ± )\ 2 



(5.13) 
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(5.14) 
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is that function of p and p± determined by the subsitutions 
(|5.6| ) and (|5.7|) . Using the relations, 



together with the definitions 
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dp 
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(5.15) 



which follow from the definition ( |5.7[ ) and Ai"(z) = zAi(z), the differentiation in ( [5.1 3| ) can be 
carried out explicitly with the result, 
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(5.16) 



where we have neglected all terms of order A -3 and higher within the curly brackets. 

The expressions (|5.13|) , (|5.14| ) and ( |5.16|) are the main results of this paper. In order to 
understand the physics of particle creation that is captured in these expressions we compare the 



lowest order asymptotic expression for the adiabatic particle number ( 5.14 ) with the analogous 
exact expression in terms of parabolic cylinder functions for constant external electric field. 
The results are plotted in Figs. 2 through 4 for A = 1, 2 and 10 respectively. We see that 
the asymptotic expansion in terms of Airy functions reproduces the behavior of the adiabatic 
particle number quite accurately, even for moderately small A of order one. The other important 
feature to notice about these figures is the relatively sharp increase in particle number right 
around u = p = 0. The transients after this particle creation event then settle down to the value 
\/3\ 2 = exp(— 7rA) which is independent of the initial longitudinal momentum. 

Thus, the exponentially small Schwinger particle creation effect is captured very well by 
the leading order term in the uniform asymptotic expansion of |/?| 2 . Notice that the uniform 
asymptotic expansion for the source term works very well even at the expected limit of its validity 
at A = 1. As a mathematical aside we remark that the exponentially small contribution to an 
adiabatic invariant quantity such as the particle number A/k has been studied by various authors 
and bounds obtained in the general case ||24j| . However, for this particular case of constant electric 
field and mode equation ( |4.3|) leading to Weber parabolic cylinder functions, it has apparently 
not been noticed that the asymptotic expansion of the solutions of this equation, uniformly 
valid on the real time axis, allows one to calculate the exponentially small secular change in the 
adiabatic invariant A/k analytically. The same observation could clearly be generalized to other 
differential equations for which uniform asymptotic expansions are known. 



24 



A=1 



0.08 



0.06 



33, 0.04 



0.02 








I- 



0.2 
0.15 : 

0.1 
0.05 : 



10 



A=1 




10 



FIG. 2. The exact (solid curve), uniform (dashed curve) and adiabatic step function (dotted 
curve) asymptotic expansions of the adiabatic particle number and its time derivative for a constant 
electric field with A = 1 and k = 0. The particle numbers approach the same value e _7r = 0.0432 as 
t — > oo, although each A4 experiences a sharp rise at a different time, viz. near zero kinetic momentum, 
p = k + eEt ~ 0. The delta function at t = in the dotted curve of the second figure obtained from 



differentiating (5.27) is not shown. 
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FIG. 3. Same as Fig. 2 but for A = 2. The particle numbers approach the same value e~ 27r = 
0.00187 as t — ► oo, although each A4 experiences a sharp rise at a different time, viz. near zero kinetic 
momentum, p = k + eEt ~ 0. 
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A=10 
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FIG. 4. Same as Figs. 2 and 3 but for A = 10. In this case the magnitude of the step at late 
times, e~ 10n = 2.27 x 1CT 14 is much smaller than the transient effects visible in the plot, and all three 
curves are very nearly (anti)symmetric around t = 0, showing that a nearly equal number of particles 
is created and destroyed. As in the previous figures the delta function at t = in the dotted curve of 
the second figure is not shown. 
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The sharpness of the creation event at u = is clearly determined by the wave equation (f4.3|) 
to be Au ~ Az or 



At ~ = r d (5.17) 

eE 

which is the time scale for the growth of a sizable fraction of the final antiparticle amplitude 
in the quantum wave function. This time scale (which is also the time scale for the classical 
acceleration by the electric field to bring a charged particle to relativistic velocities) must be 
long compared to the quantum phase coherence time r qu , in order for the creation process to be 
described by a local approximation to the nonlocal Vlasov equation ( [3.16|) , i.e. 



Td ~ Xr qu > r qu . (5.18) 

Hence the Markov limit of the Vlasov equation requires weak electric fields A ^> 1 which is 
what we have assumed in the uniform adiabatic expansion of the source term. Conversely, if 
we consider the opposite limit where the electric field is strong, so many particles are created 
so rapidly in time that the individual particle creation events cannot be distinguished one from 
another during the quantum coherence time r qu . It is clear that in this case significant wave 
amplitude coherence during the creation process can be expected and we cannot hope to approx- 
imate the effects of such copious and coherent particle creation by a Boltzmann- Vlasov source 
term local in time, which takes no account of the prior time history. Indeed in this strong field 
limit these "particles" are not particles at all in the usual sense but are more accurately to be 
thought of as coherent wave amplitudes which lie outside of any classical or semiclassical kinetic 
particle description. 

Restricting ourselves then to weak fields these coherence effects do not need to be considered 
explicitly and are built into the initial conditions of the vacuum at t — — oo once and for all. 
However in analyzing the particle creation process in a constant field and deriving ( J5.16Q we have 
also assumed that the electric field does not vary over the typical time of the variation of A/"k- 
Thus in order to use Q5.16p in situations involving a time evolving electric field we also require 
that its time scale of variation r p i be much larger than the time scale of the creation event, i.e., 

t p i > r d . (5.19) 

If this second inequality holds then it should be possible to coarsen our time resolution still 
further by not attempting to resolve the time scale r c /. On these still longer time scales it becomes 
reasonable to approximate the sharp growth of the antiparticle amplitude near u = p z = as a 
step function, provided only that we account for the integrated value of the step from — oo to 
+oo. This is what we wish to explain next. 

Let us first reiterate that the uniform asymptotic expansion in terms of Airy functions is 
indeed essential to capturing the step explicitly in Figs. 2-4 and that the Schwinger effect is 
lost completely if a naive WKB expansion in powers of 1/A is used instead. This may be seen 
explicitly by taking the large A asymptotics of the Airy functions in ( |5.12|) . To this end we note 



the Airy functions may be represented in terms of Hankel functions of the first kind, 

i 1 5-n-i 1 (I) 

z-*M'(z) = -^e~^ w^H^Hw) (5.20) 
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with the branch cut of along the negative w axis, and w ranging from +00 — ^ to —00 — ^ 
along the horizontal contour displaced by — ^ from the real axis, as u ranges from —00 to +00, 
according to Q5.6|) . Taking the large A limit is equivalent to taking the large \w\ limit of the 
Hankel functions, which depends critically on the phase of w. This phase depends in turn on 
the sign of u from ( |5.6| ). When u < 0, then |argw| < tt and we can use the standard asymptotic 
expansion of the Hankel functions, 
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for large \w\ to find 
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(5.21] 
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Since w = — A£ is linear in A this shows that the linear combination of Airy or Hankel functions 
in /3k is of order A -1 and contains no A term, as stated above. 

Thus if t — > —00 with eE fixed or if A — > 00 with k + eEt < fixed, the adiabatic particle 
number vanishes. On the other hand if u — > +00 we cannot use ( 5.21 ) directly because arg 
w — > — tt in this limit, and the condition on the phase is not satisfied. Instead we must first use 
the connection formula, 
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(5.23) 



into the proper range in order to apply Q5.21 ). Then we find 



to bring the phase of w' = e 
that the e lw terms from Hj) 1 ^ again vanish like e™|w| _1 for large \w\, but that now there remains 
in addition the opposite frequency e~ lw term which gives 
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(5.24) 



which is finite as t — >• +00 with eE fixed. As A — > 00 with k + eEt > fixed this term is 
exponentially small compared to the ordinary A^ 1 contribution. 

In this way the uniform asymptotic expansion in terms of Airy or Hankel functions which 
contains the exponentially small Schwinger particle creation becomes non-uniform in time, de- 
pending on the sign of k + eEt, if the further asymptotic expansion of these functions in terms 
of exponentials exp(±iw) is taken. Only the uniform expansion in (|5.4j ) and (|5.9|) can capture 
the particle creation event, and Figs. 2-4 show that this it does quite accurately even at the 
lowest nonvanishing order of the expansion. This exercise in asymptotic expansions as well as 
the explicit behavior in time of the adiabatic particle number in Figs. 2-4 does shows that we 
might try the simple adiabatic expansion of /3k according to ( |5.2|) , but that we must then add 
back by hand the exponentially small step (3 in the vicinity of the creation event at k + eEt ~ 0, 
i.e. 



(5.25) 



where 
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is the result obtained by substituting the lowest order of the standard adiabatic approximation 
for the mode functions ( J5.2p , rather than Olver's uniform expansion in terms of Airy functions. 
The Heaviside step function could be replaced by any smooth function with the correct limits 
at t — > ±00. The point is that if the second inequality ( 5.19Q holds, then in multiparticle 



collective quantities such as the mean current, integrations over large ranges of kinetic momenta 
are involved, and it makes little difference whether the continuous rise in each individual mode's 
particle number on the momentum scale T c \jeE is taken into account, provided only that the 
integral over all momenta accurately describes which modes have gone through the creation 
process. It is only this fact and the second inequality involving the collective time scale of the 
plasma that can justify replacing the continuous rise of Ak by a step function. 

In this admitedly rather crude approximation the function /3(p,p±) of ( 5.1 4j ) in terms of Airy 
or Hankel functions is replaced by 

\P(p,P±)\ 2 « - exp (-^j d{u) cos(26) + exp (-ttA) 6{u) (5.27) 

in terms of elementary functions. This approximation to (|5.14 ) is compared to the uniform 



asymptotic expansion in the dotted curves of Figs. 2-4, where it is observed that it works better 
than might have been expected, except for the region near the creation event u ~ where it is 
clearly inaccurate. The delta function obtained by differentiating the last two terms of ( 5.27|) is 



not shown in the second halves of Figs. 2-4. Notice that the oscillations in these figures are well 
represented by the cos(20) term in ( |5.27| ), which may be interpreted as the interference between 



the usual adiabatic phase oscillations and the relatively sudden particle creation event. Thus we 
see that for numerical purposes it is probably sufficient to use the approximate form of \P(p, P±)\ 2 



in ( |5.27| ) for all p, except those in a band of size several units of \Jp\ + rn 2 centered at the origin 
where the sharp (but continuous) rise of particle number takes place. When one is integrating 
over a region of p or t that is large compared to the time scale r c ; over which the rise in particle 
number takes place, the crude approximation of this rise by a step function and its derivative by 
a delta function may be sufficient, provided only that their coefficient is fixed by the Schwinger 
formula, as in ( 5.27]) . On the other hand, in the region of p ~ the true behavior is certainly 



not discontinuous on the scale t c i and the more accurate form (|5.14j ) in terms of Airy or Hankel 
functions should be used for moderately strong electric fields. 

We have now succeeded in our main purpose, namely to analyze the time structure of the 
quantum particle creation process in the adiabatic number basis, and to capture that particle 
creation event by means of a uniform asymptotic expansion of the exact wave functions of the 
constant electric field background, without any need to analytically continue or approximate 



the nonlocal integral in ( |3.16|) . Because of the reasoning earlier in this section we can proceed 



to identify the time rate of change of the adiabatic particle number in the lowest order of this 
uniform asymptotic expansion given by ( 5.16|) or the time derivative of ( |5.27| ) as the effective 



source term in the Vlasov equation which describes quantum particle creation in slowly varying 
electric fields, starting from vacuum initial conditions at t = — 00. 

One point that still requires some discussion is the effect of changing the initial conditions 
from vacuum at t = —00 to those at some finite time t Q . Indeed, the comparison of the effective 
source term in the Vlasov description with the mean field evolution in the next section requires 
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that the initial conditions be specified at a finite initial time to, not at — oo. This means that 
the effective source term which is given by Eqns. fl5.16|) will differ from the actual source term 



in its dependence on the initial data and the correlations (or lack of them) in the initial state. 
To the extent that a Markovian approximation to the source term is justified and dephasing is 
efficient we expect that the memory effects of the initial conditions will be washed out on the 
time scale of significant particle creation, and therefore that the initial conditions will affect only 
the transient behavior of the evolution for times close to to- This can be checked in more detail. 

To examine the transient effects of the initial conditions let us consider arbitrary initial data 
on the mode functions /k(£o) an d /k(to), subject only to the Wronskian condition fl2.8|) and 
finite initial energy density. The general solution of /k(£) in a constant electric field is a linear 
combination of /(±)k, 

/k(t) = A k (t )/( + )k(*) + S k (fo)/(* + )k(*) ' ( 5 - 28 ) 

By using the Wronskian condition on the mode functions we can solve for the coefficients, A k (£o) 
and -B k (t ) in terms of the initial conditions on the mode functions in the form, 

Ak(to) = i(/k(*b)/f + )k(to) - /k(to)/(* +)k (M) 

£ k (t ) = «(/k(to)/(+)k(*o) - /k(fo)/(+)k(*o)) • (5-29) 

A specific example of initial data with finite energy density is provided by the adiabatic vacuum 
initial conditions at t = t , i.e., 



/k(*o) = /k(to) 
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2^ k (t ) 
-iu k (t ) 
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/k(*o) • (5-30) 



The second term in the time derivative of the mode function is essential to insure finite initial 
energy density and is nonzero for finite electric field at initial time to. It means that a definite 
nonzero value of the pair correlation C k is being assumed in the initial adiabatic vacuum state. 

Our previous choice of vacuum initial conditions is recovered if we let to ~~ > ~°° so that 
Ak — > 1, B k — > and fk(t) — > /(+) k (t). Retaining t finite means that the general expression for 
/3k (t) in Q3.4Q with the mode functions given by (|5.28|) should be used so that 

fat, to) = -a k (to)/k(t)(/( +) k + iwk/(+)k) - ^k(M/k(W(*+)k + «Wk/(V)k) • (5-31) 

Since A k and _B k are also given in terms of /(+) k by Eqns. ( p.29| ) one can develop the uniform 
asymptotic expansion for this /3 k (i, t ) using ( |5.4j ) and ( |5.9|) , repeating the steps leading to Eqn. 
( |5.14j ) keeping t finite. 



The resulting rather complicated expression for the source term will depend on the electric 
field value at the initial time to- This expression would incorporate the initial data of the actual 
mean field evolution problem starting at to more accurately than the simple choice of initial 
conditions, A k = 1 and B k = which we have used in the source term (|5.16|) . A good probe 
of the effect of these transient terms is the electric current which is plotted in Fig. 5 for A = 1. 
The early oscillations observed in the exact current are the effect of the initial conditions (|5.30|) . 



However the linear growth with t at late times can be understood from the simple approximation 
to the particle creation (|5.27 ) by a step function. 
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FIG. 5. The linear growth of the electric current with time in the case of fixed constant background 
electric field E = 1 and e = 1. The three curves shown are the current of the exact mode functions, the 
uniform Airy approximation to them with initial conditions at to = according to ( |5.28| - 5!30| ), and the 
simple step function ansatz of Eqn. (5.27). 



For if we start at t = to (rather than at t = — oo) with no initial particles present, then the 
actual current integrated over all longitudinal momenta at time t is dominated by the conduction 
current j con d in (|3.25|) and becomes 



2eexp(-7rA) / — 9(k + eEt)6{-k - eEt ) -> — exp (-ttA) (t-t ) (5.32) 

J 2n uJk 7r 

in one spatial dimension at late times, which grows linearly with the elapsed time T = t — t 
since the initial vacuum state was prepared. This is precisely the slope which is observed in all 
three curves in Fig. 5 at late times. The second step function involving i is necessary because 
only modes with initially negative kinetic momentum can go through a creation event at p(t) ~ 
since p(t) = k + eEt is a monotonically increasing function of t for constant positive E. It is 
in fact present in the and of ( |5.29|) through /k(^o) an d /k(^o) which involve the same 
parabolic cylinder mode functions and the similar behavior near k + eEto ~ as observed in 
Figs. 2-4. It is this dependence on the initial data which provides just the momentum window 
in ( |5.32| ) which we need to justify the replacement of the longitudinal momentum integration dk 



by the total elapsed time eET in ( |4.15| ), and which led to Schwinger's result for the decay rate. 
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Such an understanding of the linear time divergence is possible only with a detailed description 
of the time structure of the particle creation process as given here. 

The fact that the current in a constant electric field grows linearly with time is important 
for another reason. For it shows that backreaction must eventually be taken into account, and 
that simple perturbation theory must break down at late enough times for any nonzero eE, no 
matter how small. These backreaction effects can be taken into account only by a systematic 
resummation of perturbation theory, such as the large N expansion advocated in refs. Q- ||, or 
by the solution of the Vlasov-Maxwell system of equations, valid when the inequalities of time 
scales r qu <C r c \ <C r p i hold. 

When the electric fields are very weak fields (eE <C m 2 c 3 /h), particle creation is negligible, the 
linear slope in Fig. 5 is very small and even in backreaction the electric field will hardly change 
at all with time. In this case essentially all the effects on moderate time scales will be transient 
effects and one should reatin the initial condition information in and B^. In moderately 
strong electric fields (eE ~ m 2 c 3 /h) where particle creation is significant Fig. 5 shows that the 
transient effects of the initial data become unimportant before long and one might just as well 
use the simpler expression for the source term with = 1 and B^ = 0, derived previously. This 
is equivalent to replacing the electric field value the particles feel at the actual time of creation 
by one assumed to have been constant for times long before the creation takes place. In that case 
the source term does not depend on the value of the electric field at the initial time to, which 
again is reasonable provided t v \ ^> r c ;. It is the quasistationary, Markov approximation for the 
source term in ( |5.16| ) or ( |5.27| ) A^ = 1 and B^ = that we compare to the actual backreaction 



evolution of mean field theory in the next section. 

We conclude this section by remarking on the relationship between the local source term 
( |4.9| ) or its asymptotic expansion, ( |5.16| ) and the general nonlocal form ( |3.16|) derived in Section 
II. For a constant electric field starting from vacuum initial conditions (|4.9| ) and (|3.16|) must be 
identical of course. If, following Rau [[HJ one neglects the Bose enhancement factor 1 + 2A/"k and 
changes variables from t' to Ax = 2@k(t') — 20k(t) then the integral in ( |3.16|) may be rewritten 
in the form, 

d eEp f° , sinlWx) , rno , 

-r-A/k = -r-T / dx — - , \ \ cos(Xx) (5.33 
dt Auj 2 J-oo coshV(z) 

where 

. k + eEt' r— 
u = e — ■==- = v A smh.Lf(x) (5.34) 

\eE\ 

is given implicitly as a function of x by the relation, 

smh(p(x)cosh.Lp(x) + (p(x) = x H — 5 + sinh" 1 = , (5.35) 

Pl + m 2 y^Jpl + m 2 ) 

for constant electric field. 



This is similar in form to Eqns. (25) and (26) of [15], the additional sinh (f in the numerator 
of ( |5.33| ) being due to the fact that we have treated charged scalars rather than fermions in this 
paper. Thus the form of the source term plotted in the second halves of Figs. 2-4 is qualitatively 
similar to those presented in [15] by numerical evaluation of an integral similar to ( |5.33| ). However, 
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the neglect of the quantum statistical enhancement (or Pauli blocking) factor 1 ± 2A/k in the 
integrand of ( [3.16| ) is valid only in the weak field limit A ^> 1. Since that has already been 
assumed in writing ( |5.33| ) one should then properly evaluate the integral in the same limit. As 
already remarked in Section 2 there is no straightforward method of performing an asymptotic 
expansion of this integral in real time without losing the exponentially small Schwinger effect: 
integrating the cos (Xx) term successively by parts will generate the simple adiabatic expansion 
which contains no exp (— tt\) term or step function. In the case of weak fields a" 1 = A ^> 1 
this effect is exponentially small in any case, so if one simply evaluates ( [5.33j ) or its equivalent 
for fermions numerically as in Fig. 1 of ref. |15| or Fig. 4 of this work, most of the numerical 



contribution to what is plotted is contained in the first (pure P^ db ) term of ( 5.27|) , which scales 
like 1/A 3 , and not the last term which gives rise to the exponentially small delta function source 
of fllQf . Hence multiplication by the factor exp(7rA/2) in Eqns. (24) and (25) of ref. |]15[ for weak 
fields is nugatory, while for strong fields a -1 = A < 1, the neglect of the factor 1 ± 2A/k in (|5.33|) 
or Eqn. (25) of [O is not justified. 



VI. BACKRE ACTION 



The source term we have derived in ( |5.16|) for vacuum initial conditions at to = ~°° must be 



modified to include induced creation when there are particles present in the initial state. Since 

l + 2A/" k = (l + 2A k )(l + 2|/? k | 2 ) (6.1) 

for Ak particles in the initial state, the correct modified source term in constant electric field is 

A/k = S{p,p ± ; E) = (l + 2N k )S vac {p,p ± ; E) . (6.2) 

In the backreaction problem the electric field will vary with time. Now the local Markov approx- 
imation to the nonlocal Vlasov equation ( |3.16| ) consists of using the source term ( |6 .2| ) with the 



constant E replaced by E(t) and the constant Ak by A/k(t) at the local time of interest. The 
replacement of E by E(t) is justified if the electric field is slowly varying (the quasistationary 
approximation), while the replacement of Ak by A/k(£) is justified if the electric field is weak 
(A ^> 1), since from ( |6.1| ) the difference between the two is proportional to |/3k| 2 which is of 
order A -2 and higher order than the terms we have retained in the asymptotic expansion of the 
source term. In this way the statistical factor 1 ± 2A/k has effectively been removed from the 
nonlocal integral kernel ( p,16| ) and we have obtained the final form of the source term for use in 



backreaction. 

Converting the independent variables of the number distribution function from canonical 
momenta k and t to kinetic momenta p and t by the definition, 

A/" k (0 = A/(p = k - eA; t) , (6.3) 

we obtain the local Vlasov equation, 

—N(p,p±;t)+eE(t)—Af(p,p ± -,t) = S(p,pr,E) = (l + 2Af(p,p r ,t))S vac (p,p ± ;E) , (6.4) 

for spatially homogeneous fields. Spatial dependence in the distribution function could be in- 
cluded on the left side of ( |6.4| ) in the standard manner, provided it is also slowly varying in space 
compared to cr d . Together with the Maxwell equation, 
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A(t) = 2eJ [dp] ^jV(p,p ± ; t) + 2ej [dp] + 2Af(p,p ± ; t)) S vac (p,p ± ; E) (6.5) 

(|5.16| ), and the denning relations (|4.2| ), (|5.6| ) and ( |5.7|) this constitutes the local kinetic approxi- 
mation to the mean field equations. 

In order to understand the time scale associated with the variation of the electric field and 
therefore the validity of the quasistationary approximation to the source term by that for a 
constant electric field consider first the Vlasov-Maxwell system ignoring particle creation. With 
the source term set to zero, (|6.4|) can be solved in closed form, viz., 



jV(p,p±;t)=jV(p + eA(t),p ± ;0). (6.6) 
Substituting this solution into (|6.5| ) and linearizing in A(t) gives 

5A(t)-2e 2 5A(t) [ [dp] t jV(p,p ± ; 0) = . (6.7) 
J on op 

Integrating the latter expression by parts demonstrates that the potential (and therefore also the 
electric field) will oscillate with a frequency, 



u p i 



2e>J{MM {M ^U-^\, (6. f 



which is the relativistic plasma frequency. In the nonrelativistic limit uj(p, p±) can be replaced 
by m and the integral, 2 / [dp] J\f(p,p±; 0) = n simply gives the total number density of particles 
present in the initial state at t — 0. Then we recover the familiar expression —>■ e 2 n/m for 
the classical plasma oscillation frequency. 

This classical plasma frequency may be obtained as well from a linear response analysis of 
the quantum mean field equations as follows. We perturb the vacuum solution for the mode 
functions, 



/k(t) = y — e-^ k * , ^ = Vk 2 + m 2 (6.9) 
with zero electric field by writing 

h(t) = Mt)+6f k (t) (6.10) 
and expand the equations of motion to first order in 5 fa, A and A. The linearized mode equation, 

d 2 



,-,2 

dt 2 " 



— 



5fa = 2keAfa (6.11) 



can be solved by making use of the free retarded Green's function, 



in the form, 



G R (t-t';k) = S ^t^e(t-t'), (6.12) 
o; k 



5fa(t) = 2ek f dt'G R (t-t';k)A(t')fa(t')+A k fa(t) + Bj*(t), (6.13) 
Jo 
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where A k and _B k are constants of integration and Re A^ = in order to preserve the Wronskian 
condition ( |2.8|) under the perturbation. The corresponding linearized Maxwell equation is 



A = e I [dk] Uk (1 + 2iV k ) Re(/£<y/ k ) - 2eA — - 



7,3 



2e 2 /' dt' AW) [ [dk] (1 + 2iV k ) sin[2cu k (t - t')} - e 2 A(t) [ [dk] ( ^ + 



/,•-' 




where 



= 2e 



[dk] — (1 + 2iV k ) Re(5 k e- 2 ^ k *) 

CJ k 



(6.15) 



is given by the initial perturbation away from the vacuum solution. 

The most direct method of solving a linear integral equation such as ( |6.14| ) is to make use of 
the Laplace transform, 



A( S ) 



dte- st A(t) . 



(6.16) 



After some regrouping of terms the Laplace transform of (|6.14j) may be put into the form 



A(s) 



2iV k 



eV 



4 



m * (1 + 27Vk) 



at (s 2 + 4^) 



dk \Uk) 
sA(0) + A(0) + B(s) 



(6.17) 



where the right hand side depends only upon the initial data. We notice in ( |6.17|) the presence 
of the two-particle threshold at s 2 = — 4o> k for the creation of a pair of charged particles which 
would give rise to an imaginary part and damping in the linear Maxwell equation. Since the 
particles are massive this imaginary part is zero if we find an oscillatory solution of the equation 
with s = ±iu p i and u p i 2m. Such a solution is easily found by setting the expression in curly 
brackets in Q6.17 ) to zero and neglecting the s 4 term: 



UJ 



pi 



2e 



R,N 



[dk] N h 



8_ (k^ 
dk 



where 



i + I/[dk]g(l + 2Ar k 



(6.18) 



(6.19) 



In 3 + 1 dimensions the combination in ( |6.19|) is independent of the ultraviolet cut-off and the 
renormalized value of the charge depends in general on the distribution iV k . The only requirement 
on the distribution is that lu p i in ( 6.1 8|) must be much smaller than 2m, in which limit there is no 
particle creation at all and the time independent iV k in ( |6.18 ) may be identified with the particle 
density in phase space A/"(p). 

Thus the linear response analysis of the quantum mean field theory gives exactly the same 
result for the plasma frequency, provided that the classical charge e appearing in (|6.8|) is iden- 
tified with the renormalized charge of the quantum theory according to (|6.19 ). This provides 
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a consistency check with the classical Vlasov transport description of the plasma, valid in the 
adiabatic or infrared limit of slowly varying mean fields, and identifies the proper correspondence 
limit of the classical coupling with that in the underlying quantum description. 

It is uj p i that sets the time scale of the variation of the electric field in the backreaction 
problem, i.e. t p i ~ 2tt/uj p i. Hence the local, quaistationary approximation to the source term in 
Q6.4Q requires that the three time scales obey 

UplT qu < iO p iT d < 1 . (6.20) 

If one starts the evolution with zero initial particles then the distribution function M changes from 
its classical value N{p + eA(t),p±; 0) due to the particle creation effect embodied in the source 
term. Since the second inequality requires T c i/r qu ~ A ^> 1 and the source term is exponentially 
small in A, the number of created particles n and therefore the plasma oscillation frequency 
will also be exponentially small in A. Hence the time for enough particles to be produced to 
significantly influence the electric field will be exponentially long and the second inequality in 
(|6.20| ) will also be satisfied automatically. Thus our local, quasistationary approximation scheme 
for the source term is valid a posteriori, and we would expect even the cruder approximation of 
the source term by ( |5.27| ) to be not far from correct. 

Indeed in previous work we have shown that solving a Vlasov system with a phenomenological 
source term of the form 

(1 + 2M(p, P± ; t))\eE\ ln(l + e"^) 5{p) , (6.21) 

reproduces results quantitatively similar to the mean field theory calculation of charged matter 
field coupled to a classical electric field |2j . In the present work we have shown that no logarithm 
should be present in the source term, i.e. ln(l + e _7rA ) in ( |6.21|) should be replaced by simply e~ nX , 



and that in fact, the particle creation event is continuous and can only be crudely approximated 
as a sharp step function in ( |5.27| ) with some loss of information about the true time structure of 



the event, as demonstrated in Figs. 2-4. However if the second inequality in ( |6.20| ) is valid then 
the evolution of the mean electric field on the time scale of cu^ 1 should be affected but little by 
the further approximation of the source term by a delta function. 

In order to test the validity of this approximation we present numerical results for the 
Maxwell- Vlasov system of equations ( |6.4j ) and (|6.5| ) with both the new source term ( |5.16| ) and 



the old delta function source term (but with no logarithm), and compare the results to the exact 
solution of the mean field evolution of the mode functions ( |2.9| ) coupled to the Maxwell equation, 
(|2.17| ) for scalar QED. The electric field evolution is plotted in Fig. 6 for the three cases. We 



observe that the corrected delta function source term gives qualitatively correct results, but the 
new source term ( |6.2D does a better overall job of tracking the mean field evolution, particularly 
by getting a more accurate value of the plasma frequency at late times, where the old source 
term begins to drift out of phase. The new source term (|6.2|) also drifts out of phase eventually, 
but at a much slower rate, or in other words it more accurately estimates the plasma frequency 
of the collective motion. We deliberately chose moderately large values of the coupling e = 1 and 
the initial electric field E — 1 in order to amplify the small discrepancy between the mean field 
evoluition and that of the new source term. For evolutions at e = 0.1 such as in earlier work 0, 
the discrepancy is negligible on the scale of the plot. 
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FIG. 6. The evolution of the electric field in one space dimension, according to the exact mean field 
Eqns. pl7|) , the new source term ( |6.2| ) derived in this paper, and the old source term used previously 
( |6.2l| ), but with no logarithm, for initial electric field eE = e 2 = m 2 = 1, and no particles present in the 
initial state. The new source term tracks the mean field evolution more accurately than (6.21) which 
gives a too small plasma oscillation frequency at late times. 
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FIG. 7. The evolution of the current in one space dimension, for the same initial conditions as Fig. 
6. Both the source terms for the Vlasov equation neglect the oscillations of the current on the time 



scale T qu ~ 1. 



In Fig. 7 we display the electric current ( |3.25|) for the three evolutions. The new feature 
observed here are the rapid oscillations of the quantum mean field evolution on the time scale 
r qu and their complete absence from the evolutions with the two local Vlasov source terms 
which follow the value of the current averaged over this rapid time scale. This is in accord with 
our previous discussion of the neglect of such quantum coherence effects in any local transport 
description. The particle distribution function M is plotted as a function of k for the mean 
field and Vlasov evolutions at a particular value of t in Fig. 8. We observe the same quantum 
coherence effects in the mean field evolution here in the rapid oscillations of M in momentum 
space on a scale Ak ~ 1/ cr qu as one observes as a function of time. Again these oscillations in the 
particle distribution are absent in the Vlasov evolutions. We note also the slightly negative value 
of the distribution function in the case of the new Airy source term. This is a transient effect due 
to our setting A k {t Q ) and B k (t ) to one and zero respectively. With the more accurate source 
term computed from (|5.31|) which takes account of the initial conditions this artificial negative 
region is much smaller. It also grows less and less pronounced as time progresses, and may be 
eliminated entirely by binning the distribution in momentum bins. Some small discrepancy of 
this kind is to be expected in any truncation of the unitary field theory evolution by a local 
Vlasov source term, unless that source term is always and everywhere positive, corresponding to 
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a strictly monotonic increase of total particle number and entropy, according to ( |3.23| ). It may 
be regarded as a rough estimate of the systematic error induced by the Markov approximation 
in the source term ( |6.2| ). 

In these numerical evolutions the renormalized charge of the mean field theory was chosen 
to be = 1 in order to compare to the Vlasov evolution with unit classical charge e = 1, 
according to ( |6.18| ) and (|6.19|) . In 1 + 1 dimensions where the simulations were performed the 
charge renormalization is finite and in the vacuum is given by 



m 



m 



-r,n=o 



? + 12tt 



(6.22) 



so that the finite renormalization effects for e = m = 1 are of order 1/127T ~ 0.026 or a few 
percent in the range of the simulations shown in the figures. In the extreme weak coupling 

2 

limit t- ^> 1 where the Vlasov approximation becomes more and more accurate, this finite 
renormalization effect is completely negligible. 
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FIG. 8. The particle distribution N as a function of canonical momentum k at a fixed time t = 100 
for the mean field evolution (jagged curve) and the new source term (smooth curve) derived in this 
paper, for the same initial conditions as Figs. 6 and 7. The Vlasov equation with the new source term 
is approximately a smooth average of the actual mean field evolution on spatial scales of order cr qu ~ 1 . 
The slight dip into negative values of the smooth curve becomes less and less prominent at later times. 
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VII. SUMMARY AND OUTLOOK 



Based on the Hamiltonian description of mean field theory and the existence of an adiabatic 
invariant of this evolution, we identified the (lowest order) adiabatic particle number as the most 
suitable analog for the single particle distribution function of semiclassical transport theory. 
Although not unique, this definition of particle number involves the fewest number of derivatives 
(namely zero) of the frequency u^, and hence its time rate of change is most appropriate for 
identification as the source term for the Boltzmann-Vlasov equation, which is first order in time 
derivatives. Confirming this identification, the electric current in this basis has an intuitively 
appealing and simple quasiclassical form, ( |3.25| ). Since A/"k is already adiabatic order two with 
this definition of particle number, including higher order adiabatic corrections in Ak would be 
inconsistent with the use of the source term in backreaction as well, since Maxwell's equations 
are second order in time. 

Analyzing the time dependence of the mean particle number in a constant electric field 
we derived the rate of pair creation of charged scalar particles, and clarified the time scales 
involved in the particle production phenomenon. Although formally equivalent to the quantum 
Vlasov equation ( |3.16| ) and consistent with the general projection method of Zwanzig applied 
to the density matrix in the adiabatic number basis, our approach bypasses the mathematical 
difficulties inherent in the nonlocal integral representation, and does not require explicit use of 
the projection formalism. Unlike any direct formal manipulation of the nonlocal form (|3.16| ) 
or simple WKB expansions, we used a uniform asymptotic expansion for the local source term, 
which retains the Schwinger creation effect at the lowest order of the expansion. This local source 
term is not obtained by putting to zero the phase correlations in the pair creation process, but 
rather by the assumption that the actual correlations in a time varying field can be replaced by 
those present in a constant field at at t = — oo. This can only be approximately valid when the 
electric field is very slowly varying in time, so that any actual phase correlations in the initial 
state are no longer important. 

Given the hierarchy of time scales ( OO ) we showed that a simple modification of the usual 
expansion in terms of exponential functions is nearly adequate for most analysis of the collective 
plasma effects in scalar QED. The asymptotic expansion in terms of the elementary exponential 
functions modified by the step function, which leads to the ansatz ( |5.27| ), demonstrates in a 
simple way the origin of the linear growth in time in the current which makes backreaction 
essential at late times, for any nonzero coupling no matter how weak. It also shows why taking 
the pair production source term to be proportional to 8{p) was a reasonably good proposal after 
all, although the use of the logarithm in the source term of |TT| , |T2| and subsequent references 
seems to have been due to a confusion between the rate of particle creation and the vacuum 
persistence probability. Using this ansatz in conjunction with flbMp explains the origin of the 
Bose enhancement source term, which was incorporated in the phenomenological source term 
(|6.21| ) for physical reasons. The source term obtained in explicit form from the mode functions 
in constant electric field ( p.2[ ) is in better agreement with the mean field evolution than the 
phenomenological source term ( |6.21|) , even for quite large electric fields, although the difference 
between the two is not dramatic. 

The methods employed in this paper can be readily extended to other situations of interest, 
such as fermions, or the creation of massive particles by strong gravitational fields in an early 
universe context. The limit in which such processes can be described by a semiclassical source 
term in a transport approach should be clear from the present work: one requires a clean sep- 
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aration of the three time scales r qu associated with quantum phase oscillations, t c \ associated 
with one particle creation amplitudes, and t v \ associated with the collective motion of the mean 
field(s). Conversely, it should also be clear that when such a clean separation does not exist the 
methods of this paper cannot be applied, and very likely, no semiclassical transport approach 
is appropriate or possible. Unfortunately, this includes the cases of most interest in QCD, rela- 
tivistic heavy-ion physics and early universe cosmology, where light or strictly massless degrees 
of freedom play an important role. If m = then the low momentum modes will never behave 
like classical particles admitting a Blotzmann-Vlasov description. Even pions are light enough to 
cause the hierarchy of time or momentum scales in (|6.20|) to break down in heavy-ion collisions. 
In the formation of disoriented chiral condensates the infrared instability of the low momentum 
modes and growth of a large condensate field by coherence effects is precisely the point. In cases 
such as these where Bose condensation plays a central role, the frequencies become small 
or even imaginary, the turning point (s) of the adiabatic particle number approach or reach the 
real time axis, r qu becomes large, and no simple quasuclassical particle interpretation within 
the Boltzmann framework is possible. Complementary coherent classical field methods can be 
developed in this regime, matched to a transport description of the hard modes on a case by case 
basis, but only the full field theoretic approach is powerful enough to encompass all the various 
cases in a comprehensive fashion. 
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APPENDIX A: DENSITY MATRIX IN THE ADIABATIC PARTICLE BASIS 



In this Appendix we derive the form of the Gaussian density matrix ( 2.22 ) in the adiabatic 
number basis. Since in the case of a spatially homogeneous mean electric field the density matrix 
is a product of Gaussians for each wave number k, we consider a single wave number and drop 
the subscript k in the derivation in order to simplify the notation of this Appendix. 

For each wave number we have positively and negatively charged modes obeying the time 
dependent harmonic oscillator equation (|2.9p . Because of this and using Eqns. ( |3.6| ) and ( |3 . 7| ) 
the adiabatic particle basis is that which diagonalizes the Hamiltonian of the two-dimensional 
harmonic oscillator, 



H„ 



- (Vtt + wVV + h - c ) = - ( aia + aa] + b j b + bb^ 



(Al) 



in the complex representation. The states which diagonalize this Hamiltonian are labelled by 
two quantum numbers n+ and n_ with energy u(n + + n_ + 1). In real coordinates, 



(A2) 



we can label the states by the radial quantum number n = n + + n_ and the angular quan- 
tum number m = n + — n_ corresponding to the eigenmodes of the two-dimensional harmonic 
oscillator, 



2r dr dr 



I d 2 1 \ 

H — u 2 r 2 (r6\nm) = uj(n + l)(rO\nm) 



2r 2 d6 2 



(A3) 



in polar coordinates. As is well known these wavefunctions are given in terms of the associated 
Laguerre polynomials L°(x) in the form, 







\nm) = 









f \ 2 ime ~wr 2 /2. 
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I n—m \ | 
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The normalized eigenstates themselves may be written in the form, 



\nm) 



fat 



\-2— tf)—2— 



I n+m 1 1 I n—m \ ' 
\ 2 )• \ 2 J 



t|0) 



(A5) 



with m taking on the values — n + 2k, k — 0, 1, . . . , n, so that n ± m is an even integer. 

In order to transform the density matrix from the coordinate basis to the adibatic number 
basis it is easiest first to define the coherent states, 



\sx) = exp (ia) se %x — i$ se tx ^j |0) 



n p —imx 



s e 



— \nm) 



(A6) 
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where the prime on the sum over m denotes that m is incremented by even integers. Upon 
subsituting the explicit wave functions ( |A4j) we find the wave function of these coherent states 
can be expressed in the form, 



(r6\sx) = f-)%- 2 / 2 £ ± Ve^)^-L^K 2 ). (AT) 

V7F/ n=0m=-n \~2 ) 2 

The sums in this expression may be performed in closed form by first switching the orders of the 
n and m = —n + 2k sums, and making use of the summation formula [p5|, 

oo oo 

Y, z n L-™t 2k (x) = z k Y z n L k - n (x) = z k (l + z) k e~ xz . (A8) 

n=k n=0 

The remaining sum over k from to infinity is then a pure exponential and easily performed 
with the result, 

(rO\sx) = (^) 2 exp j~7^ 2 + 2irsy/ui cos(x -9) + s 2 | , (A9) 
or in two-component vector notation, 

(f \s) = 2 exp j-"^ 2 + ZiV^r - s + s 2 | . (A10) 

This U(l) invariant exponential form may be verified also as the solution of the differential 
equation, 

(r8\H osc \s X ) = (-Y r t r lL ~ ~hw + l^) = U [ S Ta + M|sx) (AH) 



obeying the initial condition, 



i 

(,-^|.s = 0) (^)V^ 2 (A12) 



which follows from the Schrodinger equation (|A3|) and the definition of the coherent states (A"6). 

The utility of the coherent state basis is apparent from the simple exponential form of ( |A9| ) 
or ( |A10| ), since the transformation of the density matrix from the original coordinate basis, 



\ r \P\ r / = 2^2 6XP | ^ ^ ~ T ' 4£ 2 r ' r \ ^ ' 

to the coherent state basis becomes a straightforward exercise in the integration of a product of 
Gaussians, viz., 

(s'\p\s) = / d 2 f" / d 2 r(s'\f') (r'\p\r) (r\s) 



exp {^e-V + p»s> 2 + • r ^ (A14) 
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where the real coefficients A, B, C, and d are given by 



A sin •& = —2ujr]^ 



2 



B = u 2 i 4 + 1 2 V 2 + r/ 2 £ 2 + — 

C= (<r 2 -l)u;£ 2 . (A15) 

With this result in hand all that remains to be done is to expand the coherent state density 
matrix ( A14 ) in powers of s and s' to identify the matrix elements of p in the adiabatic particle 
number basis via 

oo f s m' s n e im(x'-x) 

(s'x'\p\sx)= J2 £ r (n'm\p\nm) . (A16) 



n',n=0 m ^2^±S j^S!=2 WlL±™W« 



— m \ | 
2 



The fact that the coherent state density matrix is a function of only x' ~ X an d hence only 
m! = m matrix elements of p appear in the sum is a result of the U(l) invariance of the density 
matrix. We also note that in the pure state case, a = 1, C = and the last dot product cross 
term in (|A14j) vanishes, and with it all dependence on x' — X- m that case only m = and even n 
and n' appear in the expansion, and hence the only nonvanishing matrix elements of the density 
matrix are between uncharged states with n + = n_. Conversely, if a > 1 this is no longer the 
case and p has nonvanishing matrix elements also with charged particle states with m ^ 0. 
We first expand the exponential of the dot product, 

exp^-^ = f) i m J m (-»§"»') d m{x '- x) 

oo co I \p / C \ m + 2 P 

Multiplying this by the expansion of the exponentials of s 2 , 

asp (^) = g_l_^) , e ^, (A18) 

and likewise for s' 2 yields a fourfold sum over l,l',m and p. Collecting the powers of s and s' 
by defining new summation variables, n = 21 + m + 2p and n' = 2V + m + 2p we observe that 
£ + p — ^y 11 > and V + p = - ~ m > so that m < n and m < n'. Also from the presence of 
the T function in the denominator of ( [A17p we observe that p + m must be nonnegative, which 
implies - ^ m > and - t, m > 0. Hence m > —n and m > —n' as well, and we can write 



exp i ^s 2 + +^s' 2 + ^s- s' 



Ml 



,\m+2p / s2±2_-m-2p 



n ^oA ^p!r(p + m + i)r(^- P +i)r(^- P +i)' (AiJj 
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where M = min(n, n'). We also note that pp and n '± m are necessarily integers in this expres- 



sion. 



Because of the Y functions in the denominator the final sum over p in (|A19 ) terminates at 
p = min(n, n'). However it is convenient to retain the formal infinite range of p and make use of 
the relation for the V function, 



Y{l-z) 



T{z) 



sin [TTZ 



7T 



(A20) 



for z = p- !i y 2i ,p- a +p, — — p and - IL -p i , temporarily continuing pp and to noninteger 
values to avoid the appearance of divergences in the intermediate steps. In this way the sum 
over p is recognized as the expansion for the hypergeometric function 2 F\ = F, 



■m _\p 
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and we secure 
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The finite sum represented by the hypergeometric function with integral indices may be expressed 
in terms of Jacobi polynomials Pjf^ if desired, through the relation |26[ 



f^, m ) fl + Z 
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Comparing ( |A22| ) to the general form ( |A16| ) we may identify the matrix elements of the 
density matrix in the adiabatic particle basis to be 



(n' m\p\nm) = 

~B~ VBJ \2~a) ml 



(A24) 



ri+rn \ | / n +m \ | 
I n—m \ | / n'—m \ | 

v 2 J- V 2 



„/m-n m-n' C 
F — — ,— — ;t77+1; 
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' AA 2 



This is the desired result. It may be expressed in terms of the magnitude and phase of the 
Bogoliubov transformation from the Heisenberg basis to the time dependent adiabatic particle 
basis introduced in the text. In fact, making use of the definitions ( |2.18| ), ( |3.3| ) and (|3.18| ) we 
have 
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2co>£ 2 = cx(cosh27 — sinh27 cos$) and 
2^f] = — crsinh27 sin d , 



(A25) 



and from ( |A15| ) we obtain 



A = c<j^ 2 crsinh27 and 



B = 2uj^ 2 



<Jcosh 2 7 + 



a -I- 2 



(A26) 



so that finally, 
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Since p is symmetric under charge conjugation we have 

(n' , — m|p|n, — m) = (n'mlplnm) 



(A28) 



which implies that the mean charge J2m m ( n ' m \p\ nm ) = 0- The fact that p has nonvanishing 
matrix elements with states of nonzero m implies that the fluctuations of the charge about its 
mean value is nonzero in the general case of a ^ 1. Otherwise the most important feature of the 
general result ( |A27|) for the density matrix in the adiabatic particle basis for the purposes of the 
discussion in the text is that all the off diagonal matrix elements for n 7^ n' are rapidly varying 
functions of time because of the appearance of the phase d. Since all the phase correlation 
information of the function C of ( p. 14 ) resides in these off diagonal elements, while the average 
adiabatic particle number M is sensitive only to the diagonal matrix elements of p, the Markov 
limit of quantum Vlasov equation corresponds to replacement of the density matrix by only its 
diagonal matrix elements in this basis. 

In the pure state case a — 1, F — 1, and the only nonvanishing matrix elements of p have 
m = and n = 2£, n' = 2t' both even. The general result ( |A27| ) simplies considerably in this 
case to 



(2£'m = 0\p\2£m = 0) 



sech 7(tanli7) + e 



t+t' 



(A29) 



which yields the result ( |3.19j ) quoted in the text. We should note that this expression differs from 
that used in previous work , since in the present derivation the distinguishability of the positive 
and negative charged particles was taken into account, leading to a two-dimensional harmonic 
oscillator problem with a U(l) invariance, while the expression Eqn. (15) of the second of refs. 
or Eqn. (5.24) and the entire Appendix of the third of refs. [Q was based on a single scalar 
particle species. This is appropriate for the real uncharged $ 4 theory considered in the last of 



, whereas (|A29|) is the correct expression for the charged particle case. 
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